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Abstract. Nonclassical conservation laws with viscosity arising in multiphase fluid 
and solid mechanics exhibit a rich variety of traveling-wave phenomena, including 
homoclinic (pulse-type) and periodic solutions along with the standard heteroclinic 
(shock, or front-type) solutions. Here, we investigate stability of periodic traveling 
waves within the abstract Evans function framework established by R.A. Gardner. 
Our main result is to derive a useful stability index analogous to that developed by 
Gardner and Zumbrun in the traveling-front or -pulse context, giving necessary con- 
ditions for stability with respect to initial perturbations that are periodic on the same 
period T as the traveling wave; moreover, we show that the periodic stability index 
has an interpretation analogous to that of the traveling-front or -pulse index in terms 
of well-posedness of an associated Riemann problem for an inviscid medium, now to 
be interpreted as allowing a wider class of measure- valued solutions, or, alternatively, 
in terms of existence and nonsingularity of a local "mass map" from perturbation 
mass to potential time-asymptotic T-periodic states. A closely related calculation 
yields also a complementary long-wave stability criterion necessary for stability with 
respect to periodic perturbations of arbitrarily large period NT, N — > oo. We aug- 
ment these analytical results with numerical investigations analogous to those carried 
out by Brin in the traveling-front or -pulse case, approximating the spectrum of the 
linearized operator about the wave. 

The stability index and long-wave stability criterion are explicitly evaluable in 
the same planar, Hamiltonian cases as is the index of Gardner&Zumbrun, and to- 
gether yield rigorous results of instability similar to those obtained previously for 
pulse-type solutions; this is established through a novel dichotomy asserting that the 
two criteria are in certain cases logically exclusive. In particular, we obtain results 
bearing on the nature and mechanism for formation of highly oscillatory Turing-like 
patterns observed numerically by Frid&Liu and Canic&Peters in models of multi- 
phase flow. Specifically, for the van der Waals model considered by Frid&Liu, we 
show instability of all periodic waves such that the period increases with amplitude 
in the one-parameter family of nearby periodic orbits, and in particular of large- and 
small-amplitude waves; for the standard, double-well potential, this yields instability 
of all periodic waves. Likewise, for a quadratic-flux model like that considered by 
Canic&Peters, we show instability of large-amplitude waves of the type lying near 
observed patterns, and of all small-amplitude waves; our numerical results give ev- 
idence that intermediate-amplitude waves are unstable as well. These results give 
support for an alternative mechanism for pattern formation conjectured by Azevedo, 
Marchesin, Plohr, and Zumbrun, not involving periodic waves. 
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§1. Introduction 

In this paper, we study stability of periodic traveling-wave solutions of conser- 
vation laws with viscosity 



u, f ^ M"^, B e M"^^"^, modeling flow in compressible media. Periodic solutions 
do not occur for classical systems (1.1) possessing a convex entropy in the sense 
of [Sm,Ka], for example, ideal gas dynamics or magnetohydrodynamics, being for- 
bidden by energy considerations. However, they arise in a natural way in certain 
nonclassical systems modeling media with multiple phases, for example van der 
Waals gas dynamics and elasticity, or three-phase flow in porous media, in con- 
junction with a host of other complex phenomena not seen in classical conservation 
laws: in particular, hetero- and homoclinic cycles of traveling-front and -pulse type 
solutions, and associated nonuniqueness of Riemann solutions in the corresponding 
first-order system 



see, e.g., [Sl.l-5,Sh.l-3,AMPZ.l-4,Z]. As discussed in [AMPZ.2], these features 
have the common seed of instability of constant solutions in certain regions of state 
space: in the case of phase-transitional models, regions that are "between phases." 
Such instability is frequently (but not always) associated with change in type from 
hyperbolic to elliptic in the first order system (1.2). 

In the case of van der Waals gas and solid dynamics, there is a close relation 
between the mechanical model (1.1) and the variational Calm Hilliard model for 
phase transition; see Section 7 for further details. In particular, their stationary 
theories agree, reducing to the Euler-Lagrange equations for the associated van der 
Waals/Cahn Hilliard energy: a planar Hamiltonian ordinary differential equation. 
Thus, it is no surprise that the mechanical model features the same rich solution 
structure found in the study of Cahn Hilliard equations Utt + f{u) = Uxx, or Allen 
Cahn (reaction diffusion) equations ut + f{u) — Uxx based on the same energy 
functional. Indeed, the existence of hetero-, homoclinic, and periodic cycles is 
already suggested by the planar Hamiltonian form of the stationary- wave ordinary 
differential equation. 

On the other hand, the dynamics of these models are quite different, and so a 
priori it is not clear to what extent, or in what way, this stationary structure will 
be manifested in the asymptotic behavior of solutions. In particular, one would like 
to know whether the mechanical model (1.1) indeed supports persistent transition 
layers (fronts) between phases, and, more, whether it can successfully predict, at 
least at a qualitative level, the experimentally observed phenomena of nucleation 
(spontaneous formation of transition layers in previously smooth flow) and pat- 
tern formation (e.g., Martensitic crystal structure in stressed elastic solids). In the 



(1.1) 



Ut -h f{u)x = {B{u)ux) 
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case of three-phase flow models commonly used in oil recovery simulation, with 
gas, water and oil treated as separate phases of a single fluid, both modeling and 
experimental data are less certain than in the case of the van der Waals model. 
Nonetheless, the same basic questions are of interest, now for purposes of experi- 
mental prediction and validation or invalidation of the mathematical model. These 
questions must be addressed in the context of the full, dynamical behavior of model 
(1.1), and ultimately involve the careful study of stability. 

Stability and behavior of hetero- and homoclinic, or front- and pulse-type so- 
lutions of phase-transitional models has been investigated numerically and analyt- 
ically in, e.g., [ZPM,AMPZ.l-3,LZ.l-2,GZ,Z.l-2]. For both van der Waals and 
three-phase flow models, the picture that has emerged is that traveling-front so- 
lutions connecting one pure phase to another are stable, as are simple Riemann 
patterns involving two or more traveling fronts moving away from each other with 
nonzero speed, whereas traveling-pulse solutions connecting a single phase to itself, 
and therefore Riemann patterns in which they appear, are unstable. In terms of 
time-asymptotic dynamics, the latter seem to play the role of saddle points separat- 
ing the basins of attraction of the former, attracting asymptotic states [AMPZ.l]. 
These studies essentially answer the flrst two questions posed above, showing that 
phase-transitional layers are indeed supported by the models, and can form spon- 
taneously from smooth initial data. Moreover, the experiments of [AMPZ.l] show 
that a classical, stable small-amplitude Riemann pattern not involving phase transi- 
tions may bifurcate under an initial perturbation with compact support to a stable 
large-amplitude pattern involving two or more phase-transitional layers; this is a 
stronger sense in which nucleation is seen to occur. 

The present investigation is motivated by the third question, concerning the 
possibility of pattern formation in models of form (1.1). Particularly intriguing 
are numerical experiments of Frid&Liu [FL.1-2,LF] and Canic&Peters [CP], in 
which Riemann solutions of various phase-transitional models were seen to exhibit 
highly oscillatory, Turing-like patterns reminiscent of Martensitic crystal structure 
or nucleation in phase-transitional elasticity. A natural conjecture, by analogy with 
similar phenomena in Cahn Hilliard or reaction diffusion equations, is that these are 
stable patterns consisting of fronts connecting various periodic or constant states: 
the classical mechanism for pattern formation. 

Support for this point of view was given by recent investigations of Canic sug- 
gesting a connection with Hopf bifurcation and the appearance of limit cycles [C] 
in the traveling-wave equations arising from the specific Riemann data associated 
with these phenomena. A closer look [AMPZ.4-5] reveals a rich global struc- 
ture in the phase portraits of the associated traveling-wave ordinary differential 
equations, which indeed feature large- amplitude periodic orbits resembling each of 
the observed oscillatory motifs. However, these periodic solutions appear in one- 
parameter families of varying wavelength, terminating in the infinite-wavelength 
limit at a hetero- or homoclinic cycle corresponding to a degenerate Riemann pat- 
tern of front- or pulse-type traveling-waves with common speed: a metastable. 
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slowly interacting pattern of waves which by themselves would be stable. The 
criterion that would select a particular wavelength from among such a family is 
unclear. Moreover, in further numerical experiments [AMPZ.4-5], the patterns, 
though apparently robust, do not appear to settle down into a final, periodic con- 
figuration; rather, as conjectured in [AMPZ.4], they seem to be the result of complex 
metastable interactions of an infinite pattern of slowly interacting traveling fronts 
between pure phases, driven by linear instability in the limiting states at plus and 
minus spatial infinity (the "generalized Riemann data" in the sense of [AMPZ.l]). 

More precisely, it appears that oscillations in phase originate at spatial infinity 
through a linear instability mechanism. A stationary phase calculation [AMPZ.5] 
reveals that the response under perturbation of the unstable constant states at 
infinity (i.e., the Green function of the constant-coefficient equations obtained by 
linearizing about constant solutions with those values) is well approximated by a 
sum of modulated, time-exponentially growing Gaussian wave-packets 



(1.3) 



where kj denotes the frequency for which the associated dispersion relation Xj{k) 
takes on its maximum real part (temporal growth rate), ia* := —dXj /dk{k*), 
(3* := — {1/2) cPXj / dk'^ik*), and r^, Ij denote right and left eigenvectors of the 
fiux Jacobian df{u) evaluated at the background constant state. The waves 
might be called unstable linear diffusion waves, by analogy with the picture given 
by Liu&Zeng [LZe] in the stable case. Their oscillations grow exponentially in am- 
plitude until they reach the nonlinear regime, at which time they sharpen into slowly 
moving fronts connecting approximately pure phases and do not grow further. 

Note that this process does not involve periodic waves, or even their infinite- 
wavelength heteroclinic-cycle limit. The average wavelength of the pattern is, 
rather, closely approximated by the characteristic wavelength of the linear oscil- 
lations originating at spatial infinity, i.e., the rate at which new waves are "born." 
Likewise, the "front" separating an oscillatory region from an unstable constant 
state is just the front edge of the linear Gaussian wave-packet, whose speed of 
propagation can be well- approximated by tracking the point at which the Gaussian 
envelope reaches a fixed amplitude representing the transition to nonlinear dynam- 
ics. For more detailed discussion, and a derivation of (1.3), see Appendix B. This 
description, if correct, represents a novel and nonclassical mechanism for pattern 
formation, different from that seen in Cahn Hilliard and reaction diffusion models 
for phase transition. Indeed, it seems more related to certain models for turbulence, 
in which energy from high frequency modes drives the evolution of characteristic 
large-scale structures on a lower-dimensional attracting manifold. 

It seems difficult to rigorously verify this picture of pattern formation, which by 
its nature lies outside of the usual analytical frameworks. However, it is possible 
to give indirect support by eliminating the more usual scenario involving periodic 
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waves. For, recall that the patterns of [FL.1-2,LF,CP], consist of one or more 
oscillatory regions sandwiched between two unstable constant regions extending to 
plus and minus spatial infinity. If these oscillatory regions indeed represent pieces 
of difi'erent periodic solutions, then the most likely scenario is that at least one of 
these periodic solutions should be stable; the alternative, a stable pattern consisting 
entirely of unstable pieces, would represent a new type of dynamic stabilization for 
which we know of no possible mechanism. Thus, it is strong evidence against the 
classical scenario if we can show that there exist no large-amplitude stable periodic 
solutions lying nearby the observed oscillatory patterns. 

Motivated by these considerations, we here study the general question of stability 
of periodic solutions of conservation laws with viscosity. Of particular interest is the 
situation present in the above-mentioned numerical studies, of a planar, periodic 
family of solutions, originating from a nonlinear center and bounded by a limiting 
homoclinic or heteroclinic cycle, i.e., a typical (global) Hopf bifurcation [GH,HK]. 

Our analysis is by spectral Evans-function techniques, using an analytic frame- 
work developed by R. A. Gardner [G]. Specifically, changing coordinates to a rest 
frame for the traveling wave u, we obtain the linearized equation 

(1.4) vt = Lv := {Bv:c)x - {Av)^, 

about the (now stationary) wave u — u{x). We shall investigate stability of u by the 
study of the spectrum of the linearized operator L about the wave. More precisely, 
we investigate the weak spectral stability of u{-) as a solution of (1.1), defined as 

(1.5) a{L) c {Re A < 0}, 

i.e., nonexistence of time-exponentially unstable eigenmodes. Note that failure of 

(1.5) implies exponential linearized and (appropriately defined) nonlinear instabil- 
ity; hence, the study of spectral stability is appropriate for investigation of instabil- 
ity phenomena. (Linearized stability involves different issues, and will be discussed 
in a companion paper [OZ]). 

Following Gardner [G], we study the spectrum of L "directly," via the periodic 
Evans function, a determinant D{k, A) involving the monodromy matrix of the 
linearized eigenvalue equation 

(1.6) (L -X)w = 

for L, whose zero set {k, A) corresponds to bounded solutions e^^^w{x) of (1.6), with 
w periodic; for a detailed derivation, see Section 2. This approach has the advantage 
of generality, whereas more specialized analyses based on variational [M. 1-2, LP] or 
Hamiltonian [Mi] structure may take better advantage of the specific structure of 
the underlying evolution model. For conservation laws (1.1), there does not seem 
to be any such useful structure to which we may appeal. 
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Evans function techniques have been used successfully to study the stability of 
traveling-front and -pulse solutions; however, computations in the periodic case be- 
come considerably more complicated. Indeed, to our knowledge, ours is one of only 
three such explicit computations that have been carried out using Gardner's frame- 
work, the others being low-frequency expansion of the Evans function D{k,X) in the 
wave number k, carried out by E. Eszter [Es] in a singularly perturbed limit and by 
B. Sandstede and A. Scheel [SS.l] in the large-amplitude (long-period) limit. These 
previous analyses were obtained by different techniques, in the somewhat different 
reaction diffusion setting, and under hypotheses that do not apply here: in partic- 
ular, that zero be a simple eigenvalue of the linearized operator about the wave. 
However, the main distinction of the present analysis is that it is not in principle 
restricted to any type of limiting case, giving useful stability criteria for waves of 
arbitrary amplitude or type. (On the other hand, as we shall see, such limiting 
situations can be extremely helpful in the analytic evaluation of these criteria). 
This feature distinguishes our results also from those obtained in, e.g., [Ec,Ma.l- 
2,BMi.l-2,Mi] by other than Evans function techniques, all of which concern the 
limiting cases of a large-amplitude bounding cycle or a small-amplitude constant 
solution. 

Abstract result 1. Our main result is the development of a stability index anal- 
ogous to that obtained by Gardner&Zumbrun [GZ] in the traveling-front or -pulse 
context, relating evolutionary stability to the dynamics of the traveling-wave ordi- 
nary differential equation. Specifically, we show that 

(1.7) r = sgn 7Adetc(/"(w_) > 

is necessary for stability with respect to periodic perturbations of the same period 
as the background wave, where 7 is a transversality coefficient for the traveling 
wave ordinary differential equation 

(1.8) u' = B-\u){f{u) - fiu-) - siu - «_)), 



(1.9) u{0) = u{T) = uo, 

and A := dct{dum/du-), where Um denotes the mass over one period of the periodic 
profile u{-) and U- G M"^ is an artificial parameter. For a detailed derivation of (1.8), 
(1.9) and precise definitions of 7 and A, see Sections 2 and 5, respectively. 

Similarly as in the traveling-front or -pulse case [GZ,BSZ,ZS,Z.3], the coefficient 
A is seen to be related to well-posedness of an associated inviscid Riemann problem; 
however, this must now be interpreted in the context of measure-valued solutions 
appropriate for oscillatory solutions. Alternatively, A 7^ may be viewed as ex- 
pressing the existence and nonsingularity of a local "mass map" from perturbation 
mass to potential time-asymptotic T-periodic states, an evident necessary condition 
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(by conservation of mass) for orbital stability of u with respect to nonzero mass per- 
turbations within the class of T-periodic solutions of (1.1); for further discussion, 
see Section 5. 

Abstract result 2. The stability index F detects strong (time-exponentially grow- 
ing) instabilities analogous in the front or pulse case to unstable point spectrum 
of the linearized operator about the wave. By closely related computations, we 
derive also a complementary long-wave stability criterion detecting weak (time- 
algebraically growing) instabilities, analogous in the front or pulse case to unstable 
essential spectrum, or linearized instability of the limiting, constant states con- 
nected by the profile of the traveling wave. Whereas the stability index concerns 
instability with respect to periodic perturbations with the same period as the back- 
ground wave, the long-wave stability criterion concerns instability with respect to 
periodic perturbations on a different period, in the limit as this period goes to infin- 
ity: i.e., so-called "sideband instabilities". Specifically, in the "quasi-Hamiltonian" 
case that the traveling-wave ordinary differential equation (1.8) has an integral of 
motion, we establish an illuminating small-frequency expansion 

D{k,\) = 7Adet ( - \{dum/du_)df{u_)-^ - ikT) 
+0((|/c| + |A|)"+2), 

of the Evans function, where 7 and (dum/du-) are as described just above. This 
yields a distinguished dispersion relation 

(1.11) Ao(A;) = o{k) 
and n dispersion relations 

(1.12) Xj{k) = -ikTaj + 0{k'^), j = l, ...,n, 

where aj denote the eigenvalues of {dum/du-)~^df{u-). Relations (1.12) generalize 
those seen in the constant-coefficient case, for which dum/du- reduces to /; for 
further discussion see Example 3.4 and Remarks 3.5 and 5.10. In particular, they 
show that the "generalized hyperbolicity" requirement 

(1.13) (T({dum/du-)-^df{u-)^re&\ 

is necessary for long-wave stability. 

Remark 1.1. There is a further analogy between conditions (1.7) and (1.13) 
and the stability index and long-wave stability conditions arising in the study of 
multi-dimensional planar viscous shock fronts [GZ,BSZ,ZS,Z.3], where, similarly, 
the stability index concerns strong instability with respect to planar perturbations 
respecting the symmetry of the background solution, while the long-wave stability 
condition concerns weak instability with respect to perturbations with small but 
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nonzero transverse frequency. However, the origins of the long-wave dispersion 
relations are rather different in the two cases; in particular, (1.10) appears to be 
closely tied to the assumed quasi-Hamiltonian structure, whereas the corresponding 
expansion in the shock front case is completely general. 

Numerical results. Along with these analytic tools, we develop also numerical 
techniques for approximating the spectrum of the linearized operator about the 
wave, analogous to those developed by Brin [Br.l-2,BZ] in the traveling-front or 
-pulse case. Here, the issues are rather different. For example, the evaluation of the 
Evans function in the periodic case is numerically straightforward, since it involves 
only integration of a well-behaved ordinary differential equation on a finite interval; 
by contrast, the traveling-front or -pulse case involves integration on an infinite 
interval, making the problem stiff. On the other hand, the location of the spectrum, 
that is, the zero set {k,\{k)) of the Evans function -D(-, •) becomes considerably 
more complicated than in the traveling-front or -pulse case, for which the Evans 
function depends on a single argument only. In particular, it is quite difficult to 
resolve the n -\- 1 curves (1.11), (1.12) bifurcating from the point (/c, A) = (0,0) 
Thus, our numerical techniques are effectively restricted to (A;, A) bounded away 
from the origin, and so are complementary to the analytical techniques described 
above, which include but are not limited to the Taylor expansion of D about the 
origin. 

Applications. The stability index and long-wave stability criterion are explicitly 
comparable in the same "planar Hamiltonian" case for which the index of Gard- 
ner&Zumbrun was explicitly evaluable in the traveling-front or -pulse case: namely, 
the case for which the traveling- wave ordinary differential equation (1.8) is planar 
Hamiltonian for some distinguished speed, and together yield rigorous results of 
instability similar to those obtained previously for pulse-type solutions in [GZ,Z.l] 
In this case, using a Poincare-Bendixon argument similar to that used in [GZ], we 
may express the sign of the Melnikov integral 7 appearing in the stability index in 
terms of the sign of the derivative dT/da of period T with respect to amplitude a 
in the embedding family of nearby periodic orbits of (1.8), as simply 

(1.14) sgn 7= sgn dT/da. 

This reduces both of stability conditions (1.7) and (1.13) to requirements on the de- 
rivative dum/du-, a quantity which in the small-amplitude and the large-amplitude 
homoclinic limits approaches the identity. 

In general (i.e., for intermediate-amplitude waves), dum/du^ seems difficult if 
not impossible to evaluate analytically. Nonetheless, we are able to obtain rigor- 
ous instability results through a novel dichotomy., asserting that criteria (1.7) and 
(1.13) are in certain cases logically exclusive based only on structure or symmetries 
of dum/du- rather than its precise value. The basic strategy is to show the matrix 
df{u-){dum/du-)~^ to be trace-free, in which case its eigenvalues are real if and 
only if its determinant is negative and thus F = — sgn dT/da < if dT/da > 0. 
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(Recall that df{u-){dum/du-)~^ is 2 x 2 in the planar case we consider.) Specif- 
ically, we show for the equations of van der Waals gas dynamics and elasticity 
that there exist no stable periodic solutions, of any amplitude, with the property 
that period increases with amplitude in the one-parameter family of nearby pe- 
riodic orbits: in particular that large-amplitude waves are unstable; a separate 
argument shows that small-amplitude solutions are unstable, due to Majda-Pego 
instability of the limiting constant state [MP] . In the case of a standard double- well 
potential, it can be shown that the period is everywhere monotone increasing (see 
Remark 7.2, below), and so all periodic solutions are unstable. Likewise, for the 
special (Hamiltonian) class of planar, quadratic-flux systems considered in [GZ] we 
establish instability of large-amplitude waves lying near the numerically observed 
patterns, and (by separate, continuity argument) of small-amplitude solutions. 

These two cases are quite relevant to the issue of pattern formation that originally 
motivated our investigations. The former equations, modeling phase-transitional 
fluids/solids, are precisely those that were seen to exhibit pattern formation in the 
numerical investigations of [FL.l]. The latter, which serve as qualitative models for 
multi-phase flow in porous media, are prototypical for the other models in which 
pattern formation has been observed [FL.2,CP,AMPZ.4-5]. Indeed, taken together, 
our results are strongly suggestive of nonexistence of any stable periodic waves in 
the cases where oscillatory patterns have been observed. Further support for this 
conjecture is given by our numerical experiments, in which we find instability of 
periodic solutions across their entire (numerically determined) region of existence. 
Our results thus appear to eliminate from consideration the usual paradigm of mul- 
tiple stable periodic states, pointing to a different mechanism for pattern formation 
in multiphase conservation laws than that typically seen in reaction diffusion equa- 
tions: for example, the one proposed in [AMPZ.4-5]. 

Besides the direct, physical interest of our conclusions, we point out an indirect 
contribution of this paper that is perhaps more important. Namely, we provide a 
useful, and explicit analogy between the Evans function framework for the periodic 
case and that of the better-studied traveling-front or -pulse case, different than the 
large-period limit studied by Gardner [G.l~2]. In the present paper, we have applied 
this analogy to the results of [GZ,Z.1,BSZ,ZS,Z.3] to obtain interesting instability 
results in the context of conservation laws with viscosity, i.e., necessary conditions 
for stability. A similar translation of the complementary pointwise semigroup meth- 
ods developed in [ZH] yields interesting sufficient conditions for linearized stability; 
these results will be presented in the companion paper [OZ]. 

Discussion and open problems. The quasi-Hamiltonian assumption made in de- 
riving (1.10) seems to be fundamental for the evaluation of the long-wave viscosity 
criterion. Besides appearing frequently in physical examples, quasi-Hamiltonian 
systems are shown in [OZ] to be the only type that can support periodic waves 
that are asymptotically stable in the usual, diffusive sense (see further discussion, 
beginning of Section 6). 

A related assumption on the structure of stationary solutions is made by Bridges 
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and Mielke [BMi.2] in their powerful study, by quite different, center manifold 
techniques, of instability of multidimensional but small- amplitude periodic waves 
in the Cahn Hilliard- Allen Cahn setting 

(1.15) m^u + d^t^'^xx + '^yy + dF{^), 

F, ^f, a; e ]R, y e [0,7r] with periodic boundary, m, d > 0. Namely, they assume 
that the center manifold of small bounded solutions of the stationary equation 
near some rest point is foliated by periodic solutions consisting of level sets of an 
appropriate "spatial" Hamiltonian: that is, it has the same solution structure as 
does the (planar Hamiltonian) one-dimensional stationary equation 

(1.16) = dF{^) 

obtained formally by shrinking the y-dimension to zero in (1.15). Likewise, this 
structure is present, and used in an important way in the authors' (related) land- 
mark study of Benjamin-Feir instability in [BMi.l] and in Mielke's characterization, 
using a Lyapunov-Schmidt type reduction, of linearized stability of roll-up solutions 
of the Swift- Hohenberg equation [Mi]. 

It is interesting to compare our results to those of Bridges and Mielke in [BMi.2]. 
They derive a relation (see Lemma 4.1 of the reference) 



(1.17) D{k,A) = bo{A - Ck^) + 0{\A\^ + \k\\A\ + \k\^) 
similar to our (1.10), with 

(1.18) bo^O, sgn C= sgn dT/da. 

Here, A := mX^ + dX is the generalized spectral parameter arising in the eigenvalue 



equation associated with the linearized stability problem, and D is essentially an 
Evans function for the reduced eigenvalue equation they obtain by center manifold 
reduction. From (1.17), (1.18) there follows immediately the geometric necessary 
condition 

(1.19) dT/da < 

for stability that is their main result: equivalently, a sufficient condition dT/da > 
for instability. It is to be noted that this is a sideband-type instability condition, 
and not a stability index in the sense of our condition (1.7). 

This result generalizes well-known results in the one-dimensional case (obtained 
by dropping y in (1.15)) relating dT/da > to instability, described, e.g., in [IR]. 
However, as pointed out in [BMi.2, discussion below (1.13)], dT/da < in the one- 
dimensional case also implies instability, though of a different kind: strong (large 
A) rather than sideband instability. Thus, it might well be that such instabilities 
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occur also in the multidimensional case, even though they are not detected by the 
sideband instability analysis of Bridges&Mielke. 

In this regard, it is interesting to consider the stability of periodic solutions of 
the one-dimensional version 

(1.20) m^tt + d^t = + dF{^) 

of (1.15) from the point of view of this paper. By exactly the same calculations used 
to establish (1.7), (1.10), and (1-14), we obtain for this problem a low-frequency 
expansion 

(1.21) D{k,A) = 7A - T^/c^ + 0(|Ap + |A;||A| + \k\^) 
and a stability index 

(1.22) r= sgn 7= sgn dT/da; 

indeed, the calculations reduce substantially in this case. Here, D as usual denotes 
the Evans function for the full (unreduced) eigenvalue equation associated with 
linearized stability, and A := mX^ + dX as above denotes the generalized spectral 
parameter arising in that equation. Thus, we not only recover the sideband in- 
stability analysis (1.17)-(1.19) of [BMi], but, through the stability index, detect 
also strong instabilities in case dT/da < 0. That is, we see a similar dichotomy 
in this substantially different setting to the one used here to show instability of 
solutions of (1.1); in fact, the absence of the "hyperbolic factor" A = dum/du- in 
the formulae makes the argument much simpler to apply in this case, and leads to 
a stronger result (instability without regard to sgn dT/da). Note that there is no 
requirement here on the amplitude of solutions. 

The above discussion suggests that an Evans function analysis based on "spatial 
dynamics" might be a possible direction for generalization of the multidimensional 
analysis of Bridges&Mielke, both to the case dT/da < and to the case of large- 
amplitude waves. The Evans function described by Sandstede and Scheel [SS.2], 
based on Galerkin approximation on finite subspaces, appears to be a natural can- 
didate for such investigations. This would be an extremely interesting direction for 
further study. 

The analysis of (1.20) also suggests interesting questions in the one-dimensional 
case. Our results for the vectorial van der Waals phase-transitional model are now 
seen to be natural, though somewhat weaker, generalizations of those for the scalar 
phase-transitional model (1.20). This begs the question whether one might obtain 
a weakened version of these results also for the vector version e R"^ of (1.20). 
Indeed, one obtains for more general equations muu + dut = Uxx + G{u), u e M"', 
in the quasi-Hamiltonian low-frequency expansion 

(1.23) D{k,A) =7A-5T2/c2+0(|A|2 + |A;||A| + |A;|3), 
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and a stability index F = sgn 7, where the transversahty coefficient 7 no longer 
has a simple geometric interpretation, and the coefficient S is no longer explicitly 
evaluable. Thus, we do seem to obtain a result of instability when 5 > 0, analogous 
to that obtained for the van der Waals model when dT/da > 0. However, it is not 
clear whether 6 > can in fact occur; this would be another interesting issue for 
further investigation. Likewise, in the conservation law setting, it would be interest- 
ing to determine whether or not the results for the van der Waals model carry over 
to the general planar Hamiltonian case, or even to the vector quasi-Hamiltonian 
case, without further assumptions restricting the structure of the matrix dum/du-. 
We suspect strongly that the answer is "no," but do not so far have any counterex- 
amples. 

Finally, we mention the related analysis carried out by Laugesen and Pugh [LP] 
using quite different, variational methods, of stability under periodic perturbations 
of periodic solutions (with the same period) of the thin film evolution model 

(1.24) ht = {f{h)K^ + g{h)h)^^ + ah, 

/i G M, in the zero-gravity case a = 0, which as far as we know is the only other 
treatment of stability of periodic waves in the conservation law setting (they ob- 
tain results also for a 7^ 0, but these are not relevant to our discussion). Here 
again the stationary-wave equation is a second-order nonlinear oscillator, so planar 
Hamiltonian, even though the general (nonzero speed) traveling-wave equation is 
third-order. Applying our methods to this problem, we conjecture that one should 
be able to recover a partial version of their results, namely a stability index of 
F = sgn (dT/da), where the derivative dT/da is taken with the area under h (i.e., 
the mass) held fixed, yielding instability whenever dT/da < 0.^ (Laugesen and 
Pugh in fact characterize stability by this condition, obtaining positive stability 
results as well.) What is more interesting, one might equally well hope to carry out 
a low-frequency expansion analogous to (1.10), to obtain conditions for sideband 
instability as well. These were not treated in [LP] , and do not seem to be accessible 
by their methods. 

Remark 1.2. The planar Hamiltonian structure of (1.20) is what makes possible 
the explicit evaluation in (1.21) of the coefficient of the quadratic term in k. For 
systems of general type, this can usually be done only in some asymptotic limit; 
see for example the analyses of [Es,SS.l] in the large-amplitude limit. 

Plan of the paper. In Section 2, we frame the problem, describing the equa- 
tions and assumptions under consideration, and deriving traveling-wave and eigen- 
value ordinary differential equations. In Sections 3 and 4, we define the Evans func- 
tion, following [G.1-2], and recall a result of [G.2] relating the spectra of periodic 
waves in the large-amplitude (large-period) limit to that of a bounding homoclinic; 

^The situation here is somewhat degenerate due to the fact that the traveling- wave equation 
reduces in order at speed zero; indeed, condition (H3) of section 2 is violated. However, the 
computation can still be carried out, as discussed in Remark 5.3 below. 



M. OH AND K. ZUMBRUN 



13 



we also explore briefly the more elementary small-amplitude limit, in which peri- 
odic waves approach a limiting, constant solution. In Section 5, we define a periodic 
stability index analogous to that of the traveling-front or -pulse case, and establish 
the fundamental relation (1.7). In Section 6, restricting to the quasi-Hamiltonian 
case, we derive by similar techniques the small-frequency expansion (1.10). This 
yields, in particular, an appealing formula for the "averaged," or effective constant- 
coefficient equation governing behavior under perturbation, and also gives further 
details completing the description of large- amplitude behavior given by Gardner 
in [G.2]. Finally, restricting further to the planar Hamiltonian case, and applying 
a Poincare-Bendixon argument similar to that used in [GZ], we obtain the sign of 
the Melnikov integral 7 appearing in the stability index in terms of the sign of 
the derivative dT/da of period T with respect to amplitude a in the embedding 
family of nearby periodic orbits; In Section 7, we use the analytic tools developed 
in Sections 5 and 6 to establish rigorous instability results for the two classes of ex- 
ample systems described above. In Appendix A, we describe a numerical algorithm 
for location of the spectrum, and carry out systematic numerical experiments for 
the same two classes of example system; these support and in some cases extend 
our earlier analytical results. Finally, in Appendix B, we describe in detail the 
alternative mechanism for pattern formation proposed in [AMPZ.4-5], based on 
metastable configurations of slowly interacting fronts, driven by linear instability 
in the constant states at spatial infinity. As described in the appendix, these could 
be thought of as unstable nonlinear diffusion waves generalizing the stable versions 
described in [LZe]. 

§2. Preliminaries 

Consider a conservation law (1.1) and a periodic traveling- wave solution u = 
u{x — st), of period T, satisfying the traveling- wave ordinary differential equation 



Here, and elsewhere, denotes d/dx. Integrating (2.1) from to x, we obtain a 
first-order dynamical system 



(2.1) 



{B{u)uy = f{u)' - su', 



with initial conditions 



(2.2) 



u{0) = u{T) =: wo, 
u'{Q) = u'{T) =: ui. 



(2.3) 



u' = B-\u){f{u)-su-q\ 



(2.4) 



u(0) = u{T) = uo. 
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parametrized by {q, s) e R"'"'"^, where the "total flux" 

q :=B{u)u' — f{u) + su 
=B{uo)ui - f{uo) + suo 

is a constant of motion. Notice that the map {uq, ui) — > {uq, q) is locally invertible 
so long as det B{uo) ^ 0, by the Inverse Function Theorem, hence we have lost no 
information by this reparametrization. 

To emphasize the connection with the traveling-front or -pulse case, we assume 
further that there exists some rest point U- of (2.3), i.e., q ~ f{u-) — su-, so that 
we can rewrite (2.3), (2.4) in the final form (1.8), (1.9) given in the introduction. 
We shall use this form of the equations throughout the paper. Recall that (1.8) is 
exactly the equation satisfied by a front or pulse type traveling-wave solution 

(2.5) u = u{x — st), lim u = u±, 

X—^±00 

allowing the convenient comparison of front or pulse type and periodic solutions 
within the same dynamical system framework. Note, in the planar case of our 
main interest, that such a rest point always exists within the region bounded by a 
periodic wave, so there is no loss of generality in changing to the new coordinates; 
moreover, provided that U- is a nondegenerate rest point, det{df (u-) — si) ^ 0, 
the map U- ^ q is again locally invertible. 

We make the following nondegeneracy assumptions, analogous to those made in 
[GZ,ZH] for the traveling-front or -pulse case: 
(HO) f,BeC\ 

(HI) Re a{B(u(x))) > for aU x. 
(H2) df{u-) - si invertible. 

(H3) The traveling-wave profile w is a transversal orbit of traveling-wave equa- 
tion (1.8) augmented with s' = 0, under periodic boundary conditions u{0) = u{T); 
in particular, for fixed u is, locally, the unique T-periodic solution of (1.8) up 
to translation, even allowing variation in the speed of propagation s. 

As mentioned above, (H2) concerns nondegeneracy of the parametrization (-uq, W-, s) 
of (1.8), and thus is not strictly necessary for our analysis of periodic waves below; a 
similar analysis without this assumption can be carried out in (tto, q, s) coordinates, 
with only minor expositional changes. 

Note: Assumption (H3) does not preclude the interesting case of a planar Hamil- 
tonian ordinary differential equation possessing a family of nested periodic solutions, 
since generically the orbits of the periodic family have distinct periods; this will in 
fact be the main source of our examples in Sections 6-7. 

Now, assume without loss of generality that speed s = 0, i.e., u = u{x) is 
a stationary solution of (1.1). Linearizing (1.1) about w(-), we obtain the usual 
linearized equation (1.4), where now A, B are periodic, rather than asymptotically 
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constant as in the travehng-front or -pulse context. The eigenvalue equation for L 
is likewise 



is a periodic 2n x 2n matrix. 

Through a study of the eigenvalue equations, we shall investigate the weak spec- 
tral stability of «(•) as a solution of (1.1), as defined in (1.5). That is, we shall 
investigate whether or not the linearized operator L possesses unstable eigenmodes 
A : Re A > 0. Failure of (1.5) implies exponential linear instability with respect 
to test-function initial data, as measured in any norms whatsoever; this can be 
seen by applying the evolution operator to (f)X^ fo^^ ^ sufficiently large, where (f) 
is a (merely bounded) unstable mode and {^) '■— x{^/M), with x ^ smooth 
cutoff function that is 1 on [—1, 1] and vanishes off of [—2, 2]. By the "almost-finite 
propagation speed" property of (1.4), we find that the amplification of the resulting 
solution v{t) in going from time zero to time t is of order 



as t — > oo, precluding uniform — > stability for any choice of p and q. 

Remark 2.1. It is shown in [OZ] that strong spectral stability, defined as (1.5) 

augmented with appropriate nondegeneracy conditions (conditions (D1)-(D3) of 
[OZ]), implies linearized asymptotic stability for all p > 1, with uniform 

rates of decay equal to those for the standard heat equation. 

§3. The Evans function and the Spectrum of L 

A brief calculation reveals that L has no point spectrum in L^,p < oo: Following 
[G.1-2], we introduce the monodromy matrix 



for any integer N, for any solution W of the eigenvalue equation (2.7), whence W 
can be at most bounded and not decaying at ±oo. Indeed, we have: 




^(*)IU./lb(o)|U.~t'/^-'/% 



oo 



(3.1) M(A):=*(T,A), 
where ^'(•, A) is the fundamental solution of (2.7), i.e., 

(3.2) = A(A, x)^, *(0, A) = / 
and T is the period of the coefficients. Then, 

(3.3) W{NT) = M{X)^W{0) 
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Proposition 3.1. a{L) consists, for all L^, precisely of L°° eigenvalues, i.e., X 
such that 

(3.4) det(M(A) - 7) = 0, I7I = 1. 

Proof. If M(A) has no eigenvalue of modulus 1, then there exist k "stable", 
i.e., with modulus less than 1, eigenvalues, and 2n — k "unstable," or with modulus 
greater than 1, eigenvalues, the associated normal modes decaying exponentially 
at +00 and —00, respectively, and linearly independent. We can thus construct a 
Green function Gx, as in the asymptotically constant case, to obtain a bounded 
resolvent, hence A e p{L). For further details, see [He,Z.3-4,OZ]. 

On the other hand, L°° eigenvalues can be shown to lie in any spectrum, by 
a standard limiting argument: specifically, by showing that 

||x^^</>||l./||(L-A)x^</>||i.. ^00 

as M — > cxo for any bounded (j) such that (L — — 0, where x^i^) ■= xi^/M) 
for any C°° cutoff function x that is one on [—1,1] and vanishes off of [—2, 2]. □ 

Loosely following [G.1-2], we define the Evans function 

(3.5) D{k, A) := det(M(A) - e*'^^) 

for any (/c. A) G M x C. Note that D{k,X) is clearly jointly analytic in k and A on all 
of R X C; thus, it is somewhat better behaved than the corresponding object in the 
traveling-front or -pulse case, see [GZ,ZH] for further discussion. From Proposition 
3.1, we have immediately: 

Corollary 3.2. The spectrum of L consists of the set of all A such that D{k, A) 
vanishes for some k. 

In the terminology of [G.1-2], points A satisfying (3.4) are called 7-eigenvalues; 
we will call them /c-eigenvalues, where 7 = e*'^^. The parametrization by k gives 
a more transparent analogy to the analysis by Fourier transform of the constant- 
coefficient case, see Example 3.4 just below (or, see [OZ] or [S.1-3] for a deeper 
discussion of this analogy). 

Remark 3.3. More generally, 

(3.6) M(A) = *(L)*-i(0), 

for any matrix ^ of solutions of (2.7). The coordinate-independent representation 
(3.6) is quite useful in computations. Likewise, D{-,-) is invariant under linear 
changes of coordinates. 



(3.7) ^'-^P(a;)5', 



M. OH AND K. ZUMBRUN 



17 



with P periodic. (Note: This includes changes of coordinates at the level of (2.6), 
but also more general ones at the level of phase coordinates, (2.7)). 

It is interesting to compare to the spectrum of L considered as an operator 
on periodic functions L^[0,T]. Necessarily discrete, the spectrum of L acting on 

the bounded interval [0, T] corresponds precisely to the set of {k = 0) -eigenvalues 
of L acting on the unbounded domain; conversely, {k = 0)-stability corresponds 
to stability with respect to periodic perturbations of period T. More generally, 
{k = 27rm/n)-stability corresponds to stability with respect to perturbations that 
are periodic with period nT; letting n — >^ oo, we recover stability with respect to 
general perturbations. Applying the Implicit Function Theorem to (3.4), we find 
that the spectrum of L on the real line consists typically of curves of /c-spectra, given 
by the closure of the spectra of L considered as an operator on periodic functions 
over all multiple intervals [0, nT], n an integer. 

Example 3.4. In the constant-coefficient case, Lv := Bv^x — -^Vx with A, B 
identically constant, an elementary computation yields 

D{k, A) = Il]Zi{e^'^^'^^ - e'^^), 
where //j, j = 1, . . . , 2n denote the roots of the characteristic equation 

(3.8) dei{Bi? - A// - A) = 

associated with eigenvalue equation (2.7). On the set of A such that A is diagonaliz- 
able, this is easily seen by changing coordinates to a basis in which it is diagonal, see 
Remark 3.3; recalling that this set is dense, we obtain the full result by continuity. 
Thus, the zero-set of D consists of all fc, A such that 

(3.9) = ik (mod 27Ti/T) 

for some j. Setting /j, = ik in (3.8), we obtain the dispersion relation 

(3.10) det(A - ikA - k^B) = 0, 

recovering the standard characterization of (7(L) by Fourier transform. □ 

Note that -D(0, 0) = in the example above, for any choice of A, B, since fi = is 
an n-fold root of (3.8) at A = 0; we shall see later that this holds also in the general 
variable-coefficient case. An obvious necessary condition for stability is thus 

(3.11) {d/dk)X^{0) is imaginary, 

where A*(fc) is any smooth root of D{k, A*(fc)) = bifurcating from A*(0) = 0. For, 
otherwise there would exist exponentially unstable modes A = A*(/c) for k suffi- 
ciently small (in fact for all /c > or all k < 0, depending on sgn Re {d/dk)X*{0)). 
We shall make essential use of this simple observation later on. 
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Remark 3.5. In the constant-coefficient case, relation (3.10) yields expansions 
(3.12) = — iajk + • • • , j = 1, . . . ,n, 

for the n roots bifurcating from A(0) = 0, where aj denote the eigenvalues of A. 
Thus, we obtain the necessary stability condition of hyperbolicity, cr{A) real. 

§4. A result of Gardner 

Before beginning our main analysis, we recall an interesting related result of 
Gardner [G.2] concerning a particular (rather typical) kind of large-amplitude limit 
for a family of periodic waves. Specifically, consider the situation of a family of 
periodic traveling waves {w®}, approaching as e — > to a limiting homoclinic, or 
"pulse-type," solution u^, the vertex of which is a nondegenerate (i.e., hyperbolic) 
rest point of the associated traveling-wave ordinary differential equation. Here, 
may be taken without loss of generality to be stationary (i.e., zero speed), and 
the partial differential equations they solve to vary smoothly with s, so long as 
£ ^ is a regular perturbation, with all £-equations parabolic of some fixed order. 
Included in this framework is the usual case of traveling-wave solutions of the same 
(parabolic) partial differential equation, with possibly differing speeds, approaching 
a homoclinic separatrix in the e — > limit. 

In this scenario, as e — > 0, the period of goes to infinity, while the profile on 
a single period [0, T^] approaches uniformly to an appropriate shift of the homoclinic 
profile u^, without loss of generality to u^{x — T/2). Thus, it is natural to ask 
whether the stability properties of periodic waves may be related in the large-period 
limit to those of the limiting homoclinic wave. The following result of Gardner 
shows that this is indeed correct. 

Assume as is standard [GZ,ZH] that the linearized operator L° about the homo- 
clinic wave has essential spectrum contained in Re A < 0, so that point spectrum 
determines spectral stability, (1.5). This point spectrum may be detected by the 
vanishing of a traveling-front-type Evans function 

(4.1) D{X) det(0r, • • • , <^fe , • • • , <^^), 

where {</>^, . . . , 0^} and {(|>^_^.l, ■ ■ ■ , are appropriately chosen bases of the sub- 
spaces of decaying solutions at — oo, -|-oo, respectively, of the eigenvalue equation 
(2.7) . This can be constructed in a way that is analytic in A on all of {A : Re A > 0}; 
for details, see [GZ,ZH]. 
Then, we have: 

Proposition 4.1 [G.2]. The homoclinic limit vP is weakly spectrally stable if and 
only if the family of periodic waves {u®} is "weakly spectrally stable in the limit," in 
the sense that the spectra of the associated linearized operators about satisfy 



a{L^) C {A : Re A < r]{e)}, 
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with 7/ — > as £ — > 0. In particular, if the homoclinic u is spectrally unstable, then 
so are all for e sufficiently small. 

More precisely: 

Proposition 4.2 [G.2]. Let A C C be compactly contained in {X : Re A > 0}, 

with dA contained in the resolvent set of L°. Then, for e sufficiently small, the 
number of zeroes in A of the periodic Evans function D{-, k), for any fixed k, is 
equal to the number of zeroes of the front-type Evans function D{-). 

Proposition 4.1 is an immediate consequence of Proposition 4.2 together with 
standard sectorial bounds [Pa,He,Z.3-4] restricting the spectra of elliptic operators 
in {A : Re A > 0} to a bounded domain. Proposition 4.2 is of more general 
application, applying (under suitable, mild assumptions) also to evolution equations 
that are not parabolic. From Proposition 4.2 we obtain the revealing picture that 
eigenvalues of L°, as £ increases from zero, unfold into loops of (essential) spectrum 
of L'^] looking in reverse, we see that, as £ — > 0, the radius of the corresponding 
spectral loop decreases to zero {Proof. Denoting the eigenvalue in question as A*, 
apply Proposition 4.2 on balls i?(A*, r), with e{r) — > as r — > 0.) 

In [G.2], (an equivalent version of) Proposition 4.2 was established using the 
topological index (Chern number) construction introduced in [AGJ]. Here, we 
sketch an alternative, more elementary proof. On an initial reading, the reader 
may wish to skip this proof, which is independent from and uses a different set of 
techniques than the analysis in the rest of the paper. 

Proof of Proposition 4.2. We sketch the proof in the context of (1.1), in 
which case k = n and N = 2n in (4.1), where n as in (1.1) is the dimension of 
u. Analogous computations hold in the general case. Since both periodic- and 
front-type Evans functions are analytic in A on A C {A : Re A > 0} the result 
will follow (using Rouche's Theorem) by a winding number calculation about the 
contour F = 5A, provided we establish: 

Assertion 4-3. On any compact subset F of the resolvent set of L° intersect 
{A : Re A > 0}, 

e-«-^^(-e^'=^')-I)(A;, A.) ^ C(fc, A)D(A.), 

holds as £ ^ 0, uniformly in k, A, where C{k, A) is a nonvanishing (jointly) analytic 
function, q;_(A) is analytic, and n as above is the dimension of u. 

Proof of Assertion. By the general theory of [AGJ,GZ,ZH], the basis elements 
0^ and (f)J approach exponentially as x — > -|-oo, — oo, respectively, to the stable and 
unstable subspace of the limiting coefficient matrices Aqo := A(±oo) of (2.7) (recall 
that in the homoclinic case coefficients are not periodic, but approach limiting 
values at the same rate as does the background wave u^: exponential, in the case 
assumed here that vP{±oo) is a nondegenerate rest point of the traveling- wave 
ordinary differential equation). More precisely, 

(t>i A--- A(l)- 4>^ A--- A4>- 
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and 

0++1 A • • • A </.^„ ~ ^t+l A • • • A 

where • A ■ ■ ■ A ■ denotes an n-fold exterior algebraic product, or minor, and (p^ 
denote basis elements for the manifolds of decaying solutions at ±00 of the limiting, 
constant coefficient equations 

W = A^W 

at infinity. In turn, we have 

(^i: A---AC = e«--Ff A-.-Al/-, 
^;l+iA---A0+ =e«+-F++iA---AV^+, 

where vectors V~ and respectively span the unstable and stable subspaces of 
Aoo, and a- and denote the trace of Aqo on these respective subspaces. 
Now, let the matrix of solutions in (3.6), Remark 3.3, be chosen by 

*(0) := (V;r(0),...V;n(0),V'++i(+r/2),...,V'2+ (+7^/2)), 
where 'if}^ are determined by 

^Ji-Tm, ^+(T) <^+(+T/2). 

Then, the standard estimates of [GZ,ZH] yield that det*(A,T/2) D{X) uni- 
formly in A as £ ^ 0. 

Moreover, since -D(A) was assumed not to vanish on F, we have that 

</>r(x) A • • • A <j>-{x) ^ C_(A)e"+-Vii A • • • A y+, 

as X ^ +00, which says that decaying solutions at —00 must be growing at +00, 
and likewise 



<j>+^,{x) A • • • A 0+ (x) ~ C+(A)e---V;++i A • • • A 



+ 

2n 



as a; — > — cxo; the transmission coefficients C± as in [ZH] are determined by solving 
appropriate systems of linear equations with analytic coefficients, whose respective 
determinants factor as -D(A) times a nonvanishing analytic function. Therefore, 

(4.2) i^-{T) A • • • A ^-(T) ~ C'_(A)e«+^/V++i A • • • A F+, 
and likewise 

(4.3) Vi++i(0) A ■ ■ ■ A ^^2+ (0) ~ C+(A)e-«-^/Vf A ■ ■ ■ A F". 
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Using the representation D{k, A) = det (*(T) - e*'=^*(0)) / det ^'(0) afforded by 
Remark 3.3, combined with (4.2), (4.3) and the fact that < < «_ (more 
precisely, (/)^ , . . . , decay exponentially, while • • • 02n grow exponentially), 

we therefore obtain 
(4.4) 

D{k, A) ~ det(V;r(r), . . . , ^-(T), -e''^^^^++M, • • • , -e^'^V^2+n(0))/ det *(0) 
~ C_C+(-e^'=^)'^e("— "+)^/C+e-"+^ 
= C_(-e^'=^)"e''-^^. 

On the other hand, Abel's formula and the exponential convergence A — > Aqo as 
a; — > ±oo imply that 

D{X) - det^'(T/2) 
(4.5) = CQe-("-+"+^™ det *(T) 

~ CoC_ det(F+, . . . , F+, ...,V,-), 

where 

(7o(A) := e^o°°( Tr A^(A)- Tr A(A,a;))da; ^ Q^^y 

Comparing (4.4) and (4.5), we obtain the result, with 

C{k, A) := {Co det{V+, V+, V-^„ V,-))-\ 

□ 

Remarks. 1. As noted above, the calculation in the proof of Proposition 4.2 
also shows explicitly that the /c-loops of the spectrum shrink as T — > oo to single 
eigenvalues, thus illuminating the passage from continuous to point spectrum as 
the period T goes to infinity. Though we did not do it here, the rate of shrinking 
could be quantified by a more careful version of the same argument. 

2. The final assertion of Proposition 4.1 (the one that mainly concerns us here) is 
intuitively clear from consideration of a small multiple of an unstable eigenfunction 
of the homoclinic wave given as initial perturbation of a periodic wave — this will 
grow exponentially for some time, by the property of "almost finite propagation 
speed," or approximate localization of behavior. 

3. In the more general case that the family of orbits corresponding to {w^} con- 
verges as £ — to a heteroclinic cycle, one can carry out an entirely similar analysis 
to show that the fc-spectra of the periodic operators L^, for k fixed, e sufficiently 
small, correspond approximately to the union of the spectra of the linearized op- 
erators about each of the heteroclinic raveling waves in the limiting cycle. (This 
computation is in the spirit of "multi-bump" calculations carried out for multiple 
traveling-pulse solutions in models of nerve- impulse and optical transmission). 
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Applications. In [GZ], there was derived a stability mdea; suitable for the eval- 
uation of stability of traveling-front- or -pulse-type solutions of conservation laws. 
Precisely, this index yields the parity of the number of unstable eigenvalues A such 
that Re A > of the linearized operator about the wave; evidently odd parity es- 
tablishes instability, while even parity is consistent with stability but inconclusive. 
This index was explicitly evaluated in [GZ] for the class of planar, Hamiltonian 
quadratic- flux models, and later in [Z.l] for the equations of van der Waals gas 
dynamics: precisely the models we will study here in connection with pattern for- 
mation. In both cases, homoclinic, or pulse-type, solutions were found to have odd 
index, implying instability, while heteroclinic, or front-type, solutions were found 
to have even index, consistent with stability. Though the latter result is inconclu- 
sive, there is substantial evidence from other quarters [AMPZ.l,LZ.l-2] that the 
heteroclinic waves are indeed stable. 

Combining these results with the result of Proposition 4.1, we find for these 
models that periodic solutions approaching a bounding homoclinic must be unstable 
in the large-period limit, an observation that was already made in [GZ]. This is 
suggestive, but not conclusive evidence regarding the nature of numerically observed 
pattern formation. For, though numerically observed patterns do appear to lie near 
limiting separatriccs [AMPZ.5], the distance relative to the requirements of the 
abstract theory is difficult to quantify. Also, in both cases, the patterns appear to lie 
near not only homoclinic but also heteroclinic cycles: 3- and 2-cycles, respectively. 
Presuming, as evidence suggests, that the component heteroclinic fronts of such 
a limiting cycle are stable, we find from Remark 3 above that the approaching 
periodic waves are at worst weakly unstable, in the sense that unstable eigenvalues 
must have vanishingly small real part as £ — 0. Thus, the question of stability is in 
this case much more sensitive. Indeed, our numerical experiments in Appendix A 
indicate that such periodic waves are in the large-amplitude (i.e., large-period) limit 
stable with respect to period-T/(A; = 0) perturbations, but unstable with respect 
to perturbations of some other periods/ values of k, whereas, recall, the methods of 
this section do not distinguish between different /c-values. Similar considerations 
hold in the small-amplitude limit; see Remark 4.5 below. 

These difficulties motivate our development in the following sections of a more 
direct, and general approach to stability of periodic waves. First, rather than 
using the front or pulse stability index to obtain information in the large-period 
limit, we will define an analogous stability index in the periodic case, and com- 
pute this directly, thus obtaining, at least in principle, information about small- 
and intermediate- as well as large-amplitude waves. This corresponds essentially 
to determining the Taylor expansion in A of the Evans function at the origin 
{k,X) = (0,0). Next, performing a similar but less generally applicable Taylor 
expansion in the variable k, we will obtain a complementary, "long-wave" stability 
criterion differentiating between different values of k. Together, these will turn out 
to be sufficient to obtain rigorous instability results relevant to the pattern forma- 
tion phenomena discussed above: in particular, in the case of the van der Waals 
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equations, for waves of arbitrary amplitude. 

Remark 4.4. At least formally,^ the same kind of "multi-bump" analysis as 
sketched in Remark 3, above, suggests that, for a family of homoclinic orbits {u^} 
converging to a 2-cycle limit u^, the spectra of the associated linearized opera- 
tors L'^ should converge as £ to approximately the union of the spectra of 
the component heteroclinic waves in the 2-cycle. We have conjectured that these 
heteroclinic waves are stable, in which case their point spectra consist, in the set 
{A : Re A > 0}, of single translational eigenvalues at A = 0. Pursuing this line 
of reasoning, we deduce that the approaching homoclinics should have, besides the 
obligatory translational eigenvalue at A = 0, a second eigenvalue of vanishingly 
small real part, the sign of which determines stability or instability. Thus, the 
stability of homoclinics would be quite sensitive in the limit as they approach a 
2-cycle, a result consistent with the sensitive stability we have deduced for their 
nearby periodic waves (which must with them approach to the same 2-cycle). 

Remark 4.5. It is interesting also to consider the opposite situation from that 
studied by Gardner, namely the small- amplitude limit as a family of periodic waves 
{u'^} shrinks to a single point, or nonlinear center Uc- Specifically, let us consider 
the typical situation of a Hopf bifurcation, for which the linearization 

(4.6) v' = B-^A{uc)v 

of the traveling-wave ordinary differential equation Bu' = Au about Uc possesses 
a two-dimensional center manifold corresponding to a single pair of complex eigen- 
values a± := ±ZT of the coefficient matrix B~^A{uc)- In this scenario, the period 
T'^ of u'^ converges as e — to T*^ := 27r/T, while the amplitude shrinks to zero. 

Likewise, the eigenvalue equation (2.7) converges as £ to a constant-coefficient 
equation as considered in Example 3.4, with A = A{uc) and B = B{uc)- Referring 
to (3.8), we find that, for A = 0, there is an n-fold root /i = and n remaining roots 
jji consisting of the eigenvalues of B~^A. Consulting (3.9), and noting that roots 
±ir are by definition equal to ±2ni/T^, we find, fixing /c = 0, that A = is an 
(n-|-2)-fold root of D{0, •) = 0, while all other roots are far from the origin, and may 
in fact be stable: this is the case, for example, when B = I. Thus, stability with 
respect to periodic perturbations of period T, corresponding to stability/instability 
of the roots of the restricted Evans function D(0, •), may be quite sensitive in the 

^There is a technical difficulty associated with the accumulation of the essential spectrum at the 
imaginary axis of the linearized operators about these waves; in particular, note that the argument 
of Proposition 4.1 makes essential use, specifically in the derivation of the crucial estimate (4.4), 
of the fact that cj)'^ , . . . , (j}^ are exponentially decaying for x near +00, while (f'n+i^ • • • ' ^2n 
exponentially growing for x near —00 (which is to say they decay as a; ^ —00). It is a standard 
fact (see, e.g. [He,GZ,ZH]) that this property is equivalent to the assumption that A lie uniformly 
to the right of the essential spectrum boundary of L. On the other hand, the more general bundle 
construction of [AGJ,G.l— 2] should still go through, in conjunction with analytic continuation 
into the essential spectrum of the traveling-front- or -pulse-type Evans function using the Gap 
Lemma of [GZ,KS]; see, for example, the related analysis in [DGK]. 
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small-amplitude limit. On the other hand, waves are usually quite unstable with 
respect to perturbations on different periods, corresponding to zeros of D{k, •) with 
k 0: for example, in case B = I, we have that A has imaginary eigenvalues ±zt, 
and so hyperbolicity, hence also (3.11), is violated. 

§5. The Stability Index 

Motivated by the discussion of the previous section, we now define directly a 
stability index for periodic waves analogous to the one defined in [GZ] for traveling- 
front- or -pulse-type solutions. This index, being designed to identify strong insta- 
bilities, concerns only (k = 0)-eigenvalues, i.e., instability with respect to periodic 
perturbations of the same period T as the background solution u. For this purpose 
the choice of k is expected to be somewhat arbitrary (by continuity; see also the 
large-period analysis of the previous section). And, this allows us to make a link to 
the previous analyses of heteroclinic waves, for which the Evans function has only 
the argument A. 

As in the analysis of [GZ] in the traveling-front or -pulse case, our goal is to relate 
sgn -D(0, A) for A near +cxo to sgn -D(0, A) for A near 0, for A restricted to the real 
axis, and, in turn, to relate the latter to the dynamics of the traveling- wave ordinary 
differential equation. This approach, introduced by Evans in the pioneering papers 
[E.1-4], and widely generalized in, e.g., [J,AGJ,PW], has proven to be a powerful 
tool in the stability analysis of traveling-front or -pulse type solutions. However, 
until now it does not seem to have been carried out in the periodic case. 

Large A Behavior. We begin with the large-A limit, which admits a particularly 
simple treatment. 

Lemma 5.1. As X +oo along the real axis, sgn D{0, A) {—1)^, where n as 
in (1.1) is the dimension ofu. 

Proof. By standard Garding-type (i.e., sectorial) energy estimates, L has no 
spectrum in Re A > for |A| sufficiently large. Moreover, the Evans function varies 
continuously with respect to continuous changes in the coefficients of L, by con- 
tinuous dependence with respect to initial data of solutions of ordinary differential 
equation. Thus, the quantity 

(5.1) lim sgn D{0,X), 

A— >+oo 

with A restricted to the real axis, is both well-defined and invariant under homotopy 
in L within the class of strictly elliptic operators with periodic coefficients of period 
T. Deforming L to the Laplacian L := {d/dxY via the homotopy 

eL + {i-e)L, 

6 going from to 1, we may thus evaluate (5.1) by an explicit and elementary 
computation, which we omit. This could alternatively be carried out directly, for 
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the original operator L, using a rescaling argument as in [GZ]. A related, but much 
more complicated, homotopy argument was used in [BSZ] to treat the traveling- 
front or -pulse case. □ 

Remark. Note that in this periodic context, the sign of (5.1) is an absolute 
quantity, and does not depend on any choice of coordinates. Likewise, in the 
homoclinic case, there is a natural choice of coordinates by which the sign may 
be made absolute (namely, choosing initializing bases of the stable subspace at 
+00 and the unstable subspace at —00 that together form a basis of M^"^ with the 
standard orientation; see construction of the front- or pulse-type Evans function in 
the previous section). By contrast, a significant difficulty confronted in [GZ] for the 
traveling-front case was to relate the sign at infinity to the normalization chosen at 
A = 0. Indeed, this limited the original analysis of [GZ] to the case n = 2; for the 
extension to the general case, see [BSZ,Z.3]. 

Small A Behavior. We next address the crucial small-A case. By analogy with 
the traveling-front or -pulse case, we seek to relate small-A behavior to the dynamics 
of the traveling-wave ordinary differential equation (1.8)-(1.9). 

Notice, in the present periodic context, that (1.8), (1.9) involves 2n+l param- 
eters (u-, uo, s) rather than the n + 1 parameters («_, s) of the traveling-front or 
-pulse case, since the choice of initial condition uq is completely independent of 
the critical point U-. For each choice of parameters, there is a unique solution 
^(m-,wo,s)^^-) of (1.8). We can thus define the special separation function 



Note: This is a bit different from the usual separation function in that the vanishing 
of d{-) corresponds to existence of a periodic solution of precisely period T. A 
standard Melnikov function would be based, rather, on the Poincare return map; 
here, however, we are concerned only with period T. 



(5.2) 




Variations 



(5.3) 



du/duQ ■ Cj, j = 1, ■ ■ -n, 



satisfy the linearized traveling-wave equation 



(5.4) 



Bw' = Aw, w{0) = ej. 



Without loss of generality, take coordinates such that 



(5.5) 



hence 



(5.6) 
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Likewise, the variation 

(5.7) zx := -du/ds 
satisfies 

(5.8) Bz' = Az+{u-u-), ^(0) = 0, 
while 

(5.9) Wn+j := (du/du-) ■ ej, j = l,---n, 
satisfy 

(5.10) Bw' = Aw - f{u-)ej, w(0) = 0. 
We thus have relations 

(5.11) (dd/duo) ■ Cj — [wj], 

(5.12) -dd/ds=[zi], 

(5.13) (dd/du-) ■ ej = [wn+j], 

where ej denotes the jth standard basis element in W^. 

With these definitions, we may restate (H3) in the more quantitative form: 

(H3') dd/d{uo, s) = ([ii], [wi], . . . , [wn]) is full rank, 

i.e, (recalling that [wi] = [ux] = 0) 

(5.14) 7:=det([ii],[«;2],...,K])7^0. 

That is, for a fixed period T and equilibrium u-, the orbit u{-) is locally unique up 
to translations, even allowing variation in s; moreover, it corresponds to a transverse 
intersection of the tangent manifolds at ■u(0)*^"-'"°'*^ and u(T)^'^-''^°'^^ with respect 
to variations in (uq, s). 

Then, by the Implicit Function Theorem, there is an (n + l)-dimensional surface 

(5.15) {u-,ul) ^ {u-,uo,s) 

in parameter space for which '"0''^)(.) has period T, where := uq ■ ei denotes 
the component of uq in the ei = Ux{0) direction. Fixing Uq to factor out translation 
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invariance and fix the phase, we find n directions in which connections persist. We 
can thus uniquely specify n variations in qS'^-^'^o^s).^ 

(5.16) Wn+j e Wn+j ® Span(ii,t(;2, • • -t^n), 
by the requirement 

(5.17) = 0. 

Finally, we define the "zero- viscosity stability coefficient:" 

(5.18) A := det ( j Wn+i, " ^ j ^'^n) = det{dum/du-), 
where J Wn+j denotes Wn+j{y)dy and 

(5.19) um{u-,ul):= u^^-'^°''\y)dy 

Jo 

denotes the mass over one period of the solution u('^-^'^o,s) determined by (tt_,tto) 
via map (5.15); note that this does not depend on Uq, since it is invariant under 
phase shifts. The condition A 7^ corresponds to the requirement that there exist 
no nearby periodic orbits of period T that have the same mass as does u(-), at 
least up to linear order in the perturbation of {u-,uo, s). This has an interesting 
heuristic interpretation in the spirit of [FL.1-2] as linearized well-posedness within 
a special class of measure-valued solutions, in the limit of zero viscosity, of an as- 
sociated Riemann problem having left and right states both equal to the average 
value of u over one period: namely, the class of measure- valued solutions for which 
the associated limiting sequence consists of periodic functions with a fixed ratio 
between period and viscosity. (Note that average value is preserved both under 
compact perturbations of u{-) and, by conservation of mass, under the nonlinear 
flow of (1.1).) The latter restriction comes from the fact that we are here con- 
sidering stability only within the class of periodic functions of fixed period T; its 
somewhat awkward form reflects the link via rescaling between the long-time and 
small- viscosity limits. For a careful description of measure- valued solutions and 
their relation to asymptotic behavior of (1.1), we refer the reader to [FL.1-2]. 

Alternatively, A 7^ may be viewed as the requirement that, near u, there is a 
unique periodic orbit of period T having a given mass over one period: that 
is, the "mass map" from perturbation mass to possible time-asymptotic (periodic) 
states is both well-deflned and nonsingular. Since mass per period is preserved 
under the flow of (1.1) with periodic boundary conditions, this condition is clearly 
necessary for orbital stability of u within the class of T-periodic solutions, under 
perturbations with nonzero mass. 
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In this sense, A is precisely analogous to the corresponding one-dimensional zero- 
viscosity stability coefficient A(0, 1) defined in [ZS,Z.3]^ for the traveling-front or - 
pulse case, which has the same relations to linearized well-posedness of the Riemann 
problem and nonsingularity of the mass map. With these definitions, we have the 
following fundamental relation, analogous to the one described in [GZ,ZS,BSZ,Z.3] 
for the traveling-front or -pulse case. 

Proposition 5.2. Let (H0)-(H3) hold. Then 

(5.20) D(0,A) = A"+^(-l)"det(d/(«_)-^)7A + C)(A"+2), 
or, equivalently, 

(5.21) (a/aA)'=L'(0,0) = 0, 0</c<n, 
and 

(5.22) (a/c^A)"+^L'(0,0) = (-l)'^(n+ l)!detd/(w_)-SA, 

where 7 is the transversality coefficient defined in (5.14) and A = det{dum/du-) 
is the inviscid stability coefficient defined in (5.18), (5.19). 

Remark 5.3. In case 7 = 0, but dd/d{uo,U-, s) is still of full rank, a similar 
calculation gives 

(a/aA)"+^L>(0,0) = (-l)"(n + l)!detd/(w_)-SiA, 

where 71 now denotes the determinant of an appropriately chosen alternative trans- 
verse set of vectors, and A is as in (5.18), with Wn+i, • • • , W2n chosen again to span 
the tangent manifold to the family of nearby periodic orbits, modulo translation, 
i.e., satisfying (5.17). (Of course, these cannot now be defined as in (5.16).) Thus, 
we see precisely the same relation to linearized well-posedness of Riemann problems 
as found in the case of Lax and undercompressive shock waves in [GZ,ZS,BSZ,Z.3]. 

Remark 5.4. It may happen that 7 — > but A remains bounded as {uu-,s) 
and T approach certain limiting values. In this case, we may conclude that F = for 
the limiting periodic orbit, since D, as the uniform limit of analytic functions, has 
continuous partial derivatives as well. For example, 7 = in the small-amplitude, 
constant-coefficient limit described in Remark 4.5, since in this case there is a one- 
parameter family of periodic orbits with period T = T°, but also du/du- — > / as 
this limit is approached; see Remark 5.10, below. Thus, we may conclude that F = 
in the constant-coefficient case. This is consistent with our previous observation, 

■^The relation between the stability index, linearized well-posedness of the Riemann problem, 
and the mass map was first noted in [GZ], in slightly less explicit form; see [FreZ] for related 
applications. 
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obtained by direct calculation, that D(0, ■) must vanish in this case to order (n + 2) 
and not (n + 1). 

Proof of Proposition 5.2. Defining the fundamental set of solutions VFi, . . . , W2n 
of eigenvalue ordinary differential equation (2.7) by initialization 

(8.23) (Wu....W,^m=[J,^ 

and writing Wj =: (wj, w'^Y, we obtain by Remark 3.3 the representation 



(5.24) 



3^ 

D{0, A) = det{[Wi],. . . , [W2n])/ det(M^i(0), . . . , W^2n(0)) 



With this choice of coordinates, the Wj{X) defined here agree at A = with the 
traveling- wave ordinary differential equation variations Wj defined in (5.4), (5.9), 
satisfying the (second-order) linearized traveling- wave ordinary differential equation 

(5.25) (Bw'y = (Aw)'. 
Likewise, variations zj := Wj^ are seen at A = to satisfy 

(5.26) {Bz'y = {Az)' + Wj, z{0) = ^'(0) = 0, 
and Uj := wj^^ to satisfy 

(5.27) {Byy={Ayy + 2zj, y(0) = y'(0) = 0. 
In particular, zi satisfies 

(5.28) {Bz'Y = {AzY + Ua,, z{0) = 0, 
hence 

(5.29) Bz' = Az+ (u-uo), ^(0) = 0. 
Comparing (5.29) with (5.8), (5.10), we find that 

(5.30) zi = z\ modulo span {wn-\-ii ■ ■ ■W2n)- 
Next, integrating (5.25) from to T, we find that, at A = 0, 

(5.31) S(«o)K-] = [Bw'j\ = [Awj] = AiuoYwj]; 
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T 

Wrr = win = 0. 



hence 

(5.32) [w'j] - B~^A{uq) [wj] = 
for all 1 < J < 2n. Likewise, integrating (5.26), (5.27) gives 

(5.33) [4] - B-^A{uo)[zj] = B-\uo) ^ wj, 

(5.34) [y'^ - B-'A{uQ)[yj] = 2B-\uq) I z^, 

Jo 

respectively. In the special case j = l, (5.33) gives 

(5.35) [4] - B-^A{uo)[z,] = B-\uo) f 

Jo 

Now, Taylor expanding formula (5.24) about A = 0, we obtain for small A that 
(5.36) 

D{0, A) = (-1)" det {df{u-)-^B{uo)) 

^ . VM+AN + iA2[yi] + ... , [«;2^] + A[^2n] + ---\ 

""^'[H] + X[z[] + '^X^[y[] + ■ ■ ■ , [w',J+X[z',J + ---J- 

Subtracting B~^A{uq) times the first row from the second, and using (5.33), (5.34) 

and f wi = [wi] = [w'l] = 0, we obtain 

(5.37) 

i^(0,A) = (-l)-det {df{u.)-^Biuo)) 

X[zi] H , [W2] -i , . . . , [1027!,] H \ 

X'^B-^ Jzi + --- , XB-^Jw2 + ---, XB-^ Jw2n + --- J 

= A'^+i(-l)-det(d/(«_)-Mdetf[5i]' t^^]' [^^"H 

+ e)(A^+^), 

where J denotes . Applying now appropriate column operations, we obtain from 
(5.16), (5.17), (5.30) that 

D(0,A) = A^+i(-l)"det (4f(w_)-i) 

V Hpt f [^1]' [«^2],.--, [Wn], 0, 

\} Zi, j W2,...,J Wn, J Wn+l, ... } W2r 

+ C(A^+^) 

(5.38) = A"+^(-l)"det(4f(w_)-^) 



X det 



X det([2;i], [W2],. . . , [Wn]) det ( J Wn+l, . ., J W2n) 

+ 0(A^+2) 

= A^+i(-l)"det {df{u-)-^)-fA + C»(A"+2), 
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as claimed. □ 

Remark 5.5 The solutions Wn+i, ■ ■ ■ , W2ni satisfy the inhomogeneous equations 

(5.39) Bw'^j^j = Awn+j - df{u-)ej, M)n+j(0) = 0. 

Thus, they also can be expressed in terms of the fundamental solution ip{x) of the 
n-dimensional linearized traveling-wave equation 

(5.40) Bw' = Aw 

(different from the fundamental solution ^ of the 2n-dimensional eigenvalue ordi- 
nary differential equation (2.7)) via Duhamel's principle: 

(5.41) wn+j{x) = -(^j^ iP{x)i;{y)-^dyjdf{u-)ej, 



and therefore 



(5.42) 



[Wn+j] = Wn+j{T) 

i-T 



V(T) ^{y)-^dy)df{u.)e^. 
Likewise, 

(5.43) [w^] = {^{T)-m)ej, 

(5.44) j w,^[j%{y)dy)e,, 

(5.45) [zi]^i;iT) [ ^{y)-'{u{y) - u-)dy. 

Jo 

That is, formula (5.18) for A can be interpreted as a sort of Melnikov integral, 
like the corresponding object in the homoclinic (undercompressive) case of the 
traveling-front or -pulse theory [GZ,Z.l]. 

Remark 5.6 The quantity [^i] can clearly be substituted for [ii] in (5.14), since 
the two quantities are equal, modulo Span(w^_|_i, • • ■W2n)- 

The stability index. With these preparations, we now define the stability 
index. 

(5.46) r:= sgn (a/aA)"+^D(0, 0)L»(0, +oo). 
We have immediately: 



32 



STABILITY OF PERIODIC SOLUTIONS 



Proposition 5.7. Let (H0)-(H3) hold. Then, the parity of the number of unstable 
eigenvalues of L, considered as an operator on the space of functions periodic on 
[0, T] is even if V is positive, and odd if V is negative. In particular, F > zs 
necessary for weak spectral stability as defined in (1.5). 

Proof. Evidently, -D(0, A) is invariant with respect to complex conjugation, i.e., 
-D(0, A) = l)(0. A), by the definition of This yields the familiar fact that 

the eigenvalues of the real-valued operator L, considered as acting on the periodic 
functions on [0, T], are either real or else belong to complex conjugate pairs, whence 
the parity of the number of unstable eigenvalues is equal to the parity of the number 
of unstable real eigenvalues. But, this is clearly determined by V in the manner 
stated. □ 

Remark 5.8. The strict inequality F > is necessary for strong spectral stability 
as defined in [OZ]. In particular, it is necessary for — > linearized stability for 
any p < oo, by Proposition 1.5 and condition (D3)(i) of [OZ] combined with (5.46), 
(5.10) above. This is analogous to the situation in the traveling-front or -pulse case; 
see discussion. Section 11 [ZH], of the "neutrally stable case." 

Combining the results of Lemma 5.1, Proposition 5.2, and Proposition 5.7, we 
obtain the main result of this section, an expression for F involving only geometry 
of the phase space of the traveling-wave ordinary differential equation: 

Theorem 5.9. Let (H0)-(H3) hold. Then, 

(5.47) F= sgn 7A det (m_), 

where 7 and A = det^dum/ du-) are as defined in (5.14) and (5.18), (5.19), respec- 
tively. In particular, 

(5.48) sgn 7Adet (//(«_)> 

is necessary for weak spectral stability of u. 

Remark 5.10. In the small-amplitude limit discussed in Remark 4.5, recall that 
-D(0, ■) typically has n + 2 zeroes lying near the origin A = 0, and the rest lying in 
the strictly stable complex half-plane Re A < 0. By Propositions 5.2 and 5.7, n-\-l 
of these in fact lie precisely at A = 0, with the remaining small eigenvalue being 
stable, unstable, or zero according as the stability index F is positive, negative, or 
zero. Therefore, in this case, F > is necessary and sufficient for weak spectral 
stability with respect to periodic perturbations of period T. 

Remark 5.11. With appropriate choice of the rest point it can be shown 
that the key quantity (dum/du-) in (5.47) approaches the identity in both the 
small-amplitude and the large-amplitude homoclinic limits considered earlier. For 
the small-amplitude limit described in Remark 4.5, we may choose w_ to be the 
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rest point lying near the nonlinear center Uc to which the family of periodic orbits 
converges; for the homoclinic limit, we choose U- to be the rest point lying near the 
vertex of the homoclinic. We omit the standard, but somewhat involved calculations 
giving these results. 

This has an interesting implication in the planar case, n — 2. For, in the simplest 
case that 7 is nonvanishing for the entire family of periodic waves, we find then 
by Theorem 5.9 that the sign of F depends entirely on sgn det((i/ («_)). The 
latter is negative when u- is set equal to a homoclinic vertex, since this is a saddle 
of the traveling-wave ordinary differential equation, so that det(-B~^d/(tt_)), and 
thus also det((i/(u_)) — del B det{B~^df{u-))^ is negative, as the product of a 
positive and a negative eigenvalue. But, by similar reasoning, it is positive when 
u- is set equal to a nonlinear center, since then the eigenvalues of diei{B~^df{u-)) 
are complex conjugates. Thus, combining this observation with that of Remark 
5.10 just above, we find that if large-amplitude waves are spectrally unstable with 
respect to period-T perturbations, then small-amplitude waves are necessarily stable 
in this sense. Regarding period-T perturbations, we thus have the global picture 
of a single, small real eigenvalue detaching from the origin at amplitude zero and 
initially moving into the stable complex half-plane, later crossing zero again at 
some critical intermediate amplitude to become unstable thereafter; this picture is 
borne out by our numerical investigations in Appendix A. On the other hand, we 
see that the stability index by itself is not likely to be effective across the range of 
all amplitudes, even when it is useful in the large- amplitude limit, and for small 
amplitudes is usually not useful. 

Remark 5.12. By calculations similar to those in the proof of Assertion 4.3 
(see proof of Proposition 4.2), the expression for 7 given in (5.14) may, in the large- 
amplitude homoclinic limit discussed in Section 4, be related to a corresponding 
quantity arising in the study of stability of homoclinic, or pulse-type, traveling 
waves. For n = 2, this is the quantity F defined in (3.18), Lemma 3.4 of [GZ]; a 
similar quantity arises in the general case [BSZ,Z.3]. Combining this observation 
with that of Remark 5.11 above, one can show that the sign of the periodic sta- 
bility index converges in the large-amplitude homoclinic limit, to the sign of the 
traveling-pulse-type index for the limiting homoclinic wave. We omit the associated 
calculations as lying too far from the direction of our main interest; however, see 
the related, explicit computation (6.26) in Section 6 below. 

§6. The quasi-Hamiltonian case 

We next restrict attention to the case that the traveling- wave ordinary differential 
equation (1.8) is "quasi-Hamiltonian" in the sense that it admits an integral of 
motion for each (w_,wo), with s held fixed at the base value under consideration, 
without loss of generality s = 0. This situation arises in some physical applications; 
see Section 7. Furthermore, it has a particular mathematical interest; as shown in 
[OZ] , this is essentially the only circumstance under which one can expect standard 



34 



STABILITY OF PERIODIC SOLUTIONS 



"diffusive," or asymptotic — > linearized stability of periodic solutions: in the 
generic case, linearized — > stability can hold only for p = oo, and is bounded 
at best; see [OZ], Proposition 1.5 and discussion surrounding condition (D3). 

With this additional structure, we may determine not only the low-frequency 
behavior of the Evans function with respect to A, but the entire Taylor expansion 
behavior with respect to (/c. A): 

Theorem 6.1. Let (H0)-(H3) hold, and suppose in addition that the traveling- 
wave ordinary differential equation (1.8) is quasi- Hamiltonian for the value of s 
associated with solution u. Then, for k, A sufficiently small, 

D{k, A) = 7A det ( - X{du^/du_)df{u_)-^ - ikT) 
^^■^^ +0((|A;| + |A|)"+2), 

where 7 and (dum/du-) are as defined in (5.14) and (5.18)-(5.19), respectively. 

Proof. By assumption, (1.8) admits an integral of motion H{u-,u) for each 
(w_, Wo), i-e., 

(6.2) H{u-,u{t)) = H{u-,uq) 

for all solutions u(t), where s is held fixed. Then, perturbing the periodic orbit u 
with respect to these parameters, we find that variations 

(6.3) [wj] := {d/duQ){u^''-'''°^'\T) - w("-'"°'")(0)). 



(6.4) [wn+j] := (a/aw_)(w^"-'"°'")(T) -u("-'"°'"H0)) 

must have range lying tangent to the hypersurface H{u-,u) = H{u-,uo) at the 
base values of (tt_,wo) associated with u, i.e., orthogonal to V„-ff(ti_, mq)- For, 
differentiating H{u-, u(T)) — H{u-, u{0)) = with respect to uq yields 

(6.5) = VuH{u.,uo)[du/duo], 

where, as usual, [/] denotes f{T) — /(O); likewise, differentiating with respect to 
w_ yields 

= VuH{u-,uo)[du/du-]-\-Vu_H{u-,uo)I -Vu_H{u-,uo)I 
^^'^^ =VuH{u-,uo)[du/du^]. 

Moreover, wi := ^^(O) is also orthogonal to VuH{u-,uo), since 
(6.7) = {d/dx)H{u-,u{0)) = V^if('u_, 'Uo)wx(O). 
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Noting that [W2], ■ ■ ■ , [wn] must be independent, by (H3'), we thus find that they 
are by themselves a basis for 

(6.8) Span {[wj]} © Span (wi), 

with [zi] not required. This simply reflects the fact that, by the Implicit Function 
Theorem, the full {n + l)-parameter family of nearby periodic orbits of period T 
can be generated with s held fixed; for, thanks to (6.2) only n — 1 conditions must 
be satisfied to obtain existence. 

Defining Wi, . . . , W^n again as in (5.23), but allowing k to vary, we obtain in 
place of (5.24) the representation: 

D{k, A) = det ( - dfiu^r'Biuo)) 

X det {[Wi] - {e'^^ - I)Wi{0), [W2n] - (e^'^ - /)W^2n(0)). 

Taylor expanding this about {k, A) = (0,0), we thus obtain for small k, A that 
(6.10) 

D(0, A) = det ( - df{u-)-^B{uo)) 



X 



, f [wi] + X[zi] - ikTwi{0) -\ , [w2n] + >^[z2n] - ikTw2n{0) ^ \ 

V K] + A[4] - ikTwiiO) + ■■■, KJ + A[4J - ikTw'^niO) + ■■■)■ 

Recalling (5.4), (5.9), we find that 

{w'j - B-^Awj){0) = 0, 

- B-^Awr,+,){0) - -B-\uo)df{u.)ej 

for J = 1, ... , n. Thus, subtracting B~^A{uo) times the first row from the second, 

we obtain as in (5.37) the expression 

(6.11) 

D{k,X) = det {-df{u_)-^) 

\^ X'^ J zi, Xjw2, Xjw2n + 'ikTdf{u-)en J 

+ 0((|/c| + |A|)"+2) 

= det(-d/(«_)-^)detf jt^i' [7]' j^'.li,., , ) 

^ "'^ ^ ^ \^ jZi, XjW2, Xjw2n + ikTdf{u_)enJ 

+o{i\k\ + \x\r+^), 

where in the final equality we have used (6.8) to eliminate term ikTwi{0). Applying 
column operations as in (5.38), we thus obtain 

D{k, X) = det (-df{u_)-') det (f ^2],...,^], \ 
(6.12) ^ ^ \ jy J J yj^^^ Jw2,...Jwn, M{k,X)J 

+ o{{\k\ + \x\r+^), 
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where the n x n block M is given by 

(6 13) ^""^^ ikTdf{u-)ei, J W2n + ikTdf{u-)en) 

= (dum/du-) +ikTdf{u^)I. 

Finally, factoring the determinant on the righthand side of (6.12) using its block 
triangular form yields 

D(k, A) = 7Adet (^X(dum/du-) + ikTdf{u-)^ det - df{u-)-^^ 

+o{{\k\ + \x\r+^). 

Rearranging (6.14), we obtain (6.1). □ 

Notice that result (6.1) is consistent with our previous result (3.3). As an im- 
mediate consequence, we obtain: 

Corollary 6.2. Let (H0)-(H3) hold, and suppose in addition that the traveling- 
wave ordinary differential (1.8) is quasi- Hamiltonian for the value of s associated 
with solution u. If A := det{dum/du-) ^ 0, then, for k, A sufficiently small, the 
spectrum of L consists of{n-\- 1) smooth curves: 

(6.15) Ao(A;) = o{k), 



(6.16) \j{k) = -iajkT + o{k), j = l,...,n, 
where aj denote the eigenvalues of 

(6.17) df{u.){duJdu.)-\ 

and where (dum/du-) is as defined in (5.18), (5.19). In particular, the "effective 
hyperholicity" condition 

(6.18) (j{df{u-){dum/du-)-^)ieBl 

is necessary for weak spectral stability ofu. 

If A = 0, then (6.15) holds as before. Likewise, (6.16) holds with aj := PJ^ , 

where f3j denote the eigenvalues of {dum/ du-)df {u-)~^ , for all j such that f3j ^ 0. 
Associated with vanishing j3j, however, are nonsmooth curves 

X{k) - F/'?, 
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p/q not an integer; in particular, det{dum/du-) = implies spectral instability. 

Remark. Note that F = in this (quasi-Hamiltonian) case imphes spectral, i.e., 
exponentiali-nstahility, since then det{dum/du-) = by (5.47) together with (H2) 
and (H3'). This is to be compared with the (apparently) more subtle situation of 
the general case; see Remark 6.3. 

The formulae of Corollary 6.2 may be regarded as the natural generalizations of 
the constant-coefficient formulae of Remarks 3.5 and 4.5. In particular, compar- 
ing (6.16), (6.17) to (3.12), we see that df{u-){dum/du-)~^ plays in the quasi- 
Hamiltonian case the role played in the constant-coefficient case u = uq = U- by 
A := df{u-). That is, df{u-){dum/du-)~^ may be regarded as an "averaged" or 
"effective" convection matrix for the variable-coefficient case (indeed, this is shown 
in [OZ] to hold true in a very strong sense). Note, further, that df{u-){dUm/du-)~^ 
converges in the small-amplitude limit to df{u-), since, as pointed out in Remark 
5.11, (dum/du-) — >• /. 

Likewise, in the large-amplitude, homoclinic limit considered by Gardner, we 
find again (Remark 5.11) that df{u-){dum/du-)~^ converges to df{u-), where 
U- now denotes the vertex of the homoclinic orbit, i.e., the limiting state of the 
homoclinic wave as x — > ±oo. Thus, the spectral curves described in (6.16), (6.17) 
correspond to first order in k with the spectral curves of this limiting, constant state, 
which are in turn curves of essential spectrum for the linearized operator about 
the homoclinic wave: for a detailed discussion, see the introduction of [OZ]. The 
remaining curve (6.15) corresponds to the translational (continuous) eigenfunction 
at A = associated with the spatial derivative of the profile under consideration. 

This extends to the small-frequency regime the picture described by Gardner 
(Section 4) in the large- amplitude homoclinic limit. It can be shown by an argument 
similar to that of Section 4 (together with the observations made in the introduction 
of [OZ] regarding structure of the resolvent set of the limiting homoclinic profile) 
that the former curves in fact globally approach the essential spectrum curves of 
the limiting homoclinic wave, in the sense that they approach on a ball of radius 
going to infinity with the period; thus, they extend arbitrarily far as the homoclinic 
limit is approached. We conjecture without proof that the latter curve belongs to 
a closed loop of spectra, shrinking in the large amplitude limit to the eigenvalue 
A = associated with the translational eigenfunction of the limiting, homoclinic 
wave. 

Remark. The Taylor expansion of D can likewise be carried out in the general 
(non-quasi-Hamiltonian) case, but does not seem amenable to any such simple 
interpretation. 

The planar Hamiltonian case. Finally, we consider the planar case, n = 2, 
under the assumption that the traveling-wave ordinary differential equation (1.8) is 
in fact Hamiltonian for all w_ , with s held fixed at the base speed associated with 
u, i.e.. 



(6.19) = Tr d^{B~\u)ifiu) - f{u-) - s{u - U-))). 
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This is the class from which we will take our example systems in the following 
section. Note that, for B = constant, condition (6.19) reduces to 

(6.20) = Tr {B-\df{u) - si)), 

which is evidently independent of with this observation, it is straightforward 
to construct examples. 

From (6.19), we find that the traveling- wave equation may be expressed as a 
Hamiltonian system 

(6.21) u' = B-\u){f{u) - f{u.) - s{u - u.)) = V^H{u, ti_), 

for an appropriate Hamiltonian H{u,U-), whence H{u{t),U-) is an integral of 
motion; in particular, we see that the planar Hamiltonian systems are a subclass 
of the quasi-Hamiltonian systems defined above. 

Existence of an integral of motion in the planar case implies considerable struc- 
ture of the phase portrait of the traveling-wave ordinary differential equation. For 
example, in the situation wc consider, that the a phase portrait contain a hetero- 
clinic poly-cycle, we find that the interior of the cycle must be entirely filled with 
cycles, either periodic or heteroclinic. In the simplest case that there is only a 
single nonlinear center Uc enclosed, the interior must be made up entirely of peri- 
odic orbits; see, for example, the depiction of the homoclinic case in figure 1. In 
this situation, we may globally parametrize these orbits by amplitude o, defined as 
distance from center Uc along some (fixed) curve of steepest descent if Uc is a local 
maximum, or ascent if Uc is a minimum. In more general situations, we may still 
define a local parametrization in this way, in the vicinity of any fixed periodic orbit 
u. 



0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 
-0.5 



Figure 1. The planar Hamiltonian case: bounding homoclinic cycle. 

Let us consider the period T as a function of amplitude, T = T{a). Then, we 
have: 
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Lemma 6.4. Let (H0)-(H3) hold, and suppose in addition that the traveling-wave 
ordinary differential equation (1.8) is planar Hamiltonian for the value of s associ- 
ated with solution u, with B = constant. Then, 

(6.22) sgn 7 = sgn (dT/da) 

for ^ := det([?i], [W2]) as in (5.14). 

Proof. Let Wi, W2 denote the variations defined in (5.4), (5.9), as usual normal- 
izing Wi = Ux{0) = ^1- We may locally parametrize amplitude by location of Uq 
along a curve tangent to the outward normal r] — ±62 at u{0) to the orbit {u{x)}, 
so that 

da/duQ = T] 

at 14(0), where a(tto) is the amplitude of the orbit corresponding through ttg, with 
(tt_ , s) held fixed. Differentiating the identity 



u 



{u-,uo,s) 



(T(a(uo)))-«^''-'''°''Ho) = 



with respect to uq in direction 62, we thus obtain 

[W2] + Ux{T){dT / da){ri ■ 62) = 0, 



or 



(6.23) K] = Cux (0) = Cux (T) = Get , 
with 

(6.24) sgn (C) = sgn {dT/ da) sgn det{r],Ua:{0)). 
The quantity 

(6.25) 7 = det([^i], [W2]) = Cdet{zi{T),Ux{T)) 

can therefore be conveniently calculated, using inhomogeneous Abel's equation/Duhamel 
principal as a Melnikov integral: 

rT 

7 = C / '^^det{B-\u-u^),Ux)dx 
Jo 



(6.26) 



C / det(i? ^{u — U-),Ux)dx, 
Jo 



where in the second equality we have used the Hamiltonian property Tr (B~^A) = 
constant. More precisely, this is found by setting ^{t) := C det{zi{t),Ux{t)) and 
observing that 7 satisfies 

7' = Tr {B-^A)'y + Cdet{B-\u-U-),Ux), 7(0) = 
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by (5.4) and (5.7), (5.8); for similar calculations in the traveling-front or -pulse 
case, see [GZ]. 

But, (6.26), as in [GZ], may be viewed as a contour integral and explicitly eval- 
uated using the Gauss-Green formula as 



where is the region bounded by the orbit u{-), and i] denotes the outward unit 
normal. Recalling that Tr > by (HI), and combining with (6.25), we are 
done. □ 

Remarks. 1. Formulae (6.26), (6.27) clearly show that 



as asserted in Remark 5.6, since substitution of uq for u- leaves the result un- 
changed. 

2. Comparison of formula (6.26) to that given for F in (3.18), Lemma 3.4 of [GZ] 
shows explicitly the convergence described in Remark 5.12, since the two formulae 
are formally identical, differing only in the choice of profile u. 

3. Near a bounding poly-cyclc, it is clear that dT /da > 0, indeed dT/da +oo 
as the boundary is approached. Thus, 7 > 0, and in combination with the obser- 
vation of Remark 5.11 that (dum/du-) I and det df{u-) < in the homoclinic 
case, we find that sgn F < 0, recovering the result of instability obtained by 
Gardner's technique together with the homoclinic instability results of [GZ,Z.l]. 

It is interesting to note that our periodic instability result for the homoclinic 
limit makes no requirement on the shape of the limiting homoclinic, whereas the 
general homoclinic instability results of [GZ] require that the orbit be convex in 
the vicinity of its vertex. Indeed, in the nonconvex case, the stability index of the 
homoclinic orbit is positive, consistent with stability, an apparent contradiction 
with the Remark 2 just above. In the quadratic-fiux case considered in [GZ], it was 
shown that all homoclinic orbits are (globally) convex, hence this situation does 
not arise. However, as pointed out in [Z.4], it can certainly arise in the van der 
Waals models studied in [Z.l] (described in the following section, just below). ^ 

^The possibility of nonconvex orbits was mistakenly not considered in [Z.l]; however, it can 
be seen that they are unstable despite their positive stability index, by the final remark of that 
paper. Specifically, in the notation of the reference, the zero eigenfunction Ux of the Sturm- 
Liouville operator M in the case of a locally nonconvex orbit has two nodes, so must be at least 
the third eigenfunction of M. It follows that there exists some combination v of the first two 
eigenfunctions such that J v — 0; hence v — DV for some V E L'^, and also {DV,MDV) < 0, 
giving the result. 



(6.27) 



7 = C sgn 



= C sgn 




C sgn 



det(77,w,)|0|Tr (S'^), 



det([2;i], [^2]) = det{[zi],W2]), 
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A closer inspection quickly resolves this apparent paradox. For, it is readily- 
seen that a planar homoclinic orbit that is nonconvex at its vertex saddle must 
enclose a second^ interior poly-cycle, containing the interior branches of its stable 
and unstable manifolds. Periodic orbits are thus bounded by both the interior 
poly-cycle and the exterior homoclinic, so that the argument for Remark 2 does 
not apply. Indeed, considering the simplest case that the interior poly-cycle is a 
second, locally convex homoclinic, and applying the "multi-bump" heuristics of 
Remark 4.4 and Remark 3 above it, we obtain, formally, the correct prediction of 
an odd multiplicity of unstable eigenvalues for the periodic orbit, as the sum of even 
and odd multiplicities. 

Combining the results of Theorem 5.9, Corollary 6.2, and Lemma 6.4, we obtain, 
finally, the composite stability test: 

Corollary 6.5. Let (H0)-(H3) hold, and suppose also that the traveling-wave ordi- 
nary differential equation (1.8) is planar Hamiltonian for the value of s associated 
with solution u, with B = constant. Then, provided that (dum/du-) is invertible, 
both 

(6.28) sgn (dT/da) det {df{u^){dum/du_)-^) >0 
and 

(6.29) a{df{u_){dum/du_)-^)veal 

are necessary for weak spectral stability of u, where dT j da is the rate of change of 
period with respect to amplitude, with amplitude measured in the direction of the 
outward normal to the orbit {u{x)}. If det{dum/du-) = 0, on the other hand, then 
u is spectrally unstable. (Recall, dT/da ^ by (H3') combined with Lemma 6.4-) 



§7. Calculations for example systems 

Using the complementary stability conditions of Corollary 6.5, we now derive in- 
stability results for the two classes of model, planar Hamiltonian systems that were 
considered in [Z.l] and [GZ]: van der Waals gas dynamics with artificial dispersion- 
viscosity, and the class of planar Hamiltonian models with quadratic flux functions. 
These are directly relevant to the issue of oscillatory pattern formation, being pro- 
totypes for the two kinds of systems in which these patterns have been observed 
numerically. 

Van der Waals gas dynamics. The viscous-capillary p-system. 



(7.1) 



vt-Ux = eiVxx 
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with nonmonotone stress relation has been studied by several authors as a model 
for dynamical phase transitions in compressible van der Waals fluids/solids under- 
going isothermal motion (sec, for example, [Ja,Sh.l-4,Sl.l-5,ST]). Here, £i > 
and £2 > are related to the coefficients of viscosity and capillarity of the medium 
(£1 + £2 and £i£2, respectively), v denotes specific volume/strain, u :— u -\- SiVx 
is the velocity of the medium, and p{v) denotes pressure/stress. A typical stress 
relation is p = —W'{v)^ where 

(7.2) W{v) = \{l-v^f 

is the standard "double-well" potential. 

In physical, {v^u) coordinates, equations (7.1) take the form 

Vt-Ux = 0, 

(7.3) 

Ut + P{v)x = {Sl + £2)Uxx - {ei£2)Vxxx-, 

which is sometimes more natural for computations. 



Existence. The traveling- wave ordinary differential equation associated with 
(7.1) (for periodic and front- or pulse-type waves alike) is 



(7.4) 



^ 

£1 

£2 



—U + U- 



— s 



V — V- 
U — U- 



where s denotes the speed of propagation of the wave. For speed s = 0, this becomes 
a Hamiltonian system with Hamiltonian 



(7.5) 

where 
(7.6) 



1 1 

H{vj w, f_, U-) := {u — U-Y -P(^) 

2£i £2 

= -^{u-u-f - —{W{v-)-W{v) + W'{v-){v-V-)), 
2ei £2 



P{v) :— / {p{v-) — p{z)^ dz. 

Jv- 



For nonmonotone p, the orbits of (7.4), corresponding to level sets of H, will 
on some range of f _ include one-parameter families of periodic orbits bounded by 
heteroclinic or homoclinic cycles. For example, in the case of (7.2), there appears for 
= ±1 a 2-cycle of heteroclinic orbits connecting {v-,U-) to {—V-,U-) in either 
direction, respecting the vertical symmetry of H about the line u = U-] within this 
cycle is a one-parameter family of periodic orbits converging in the small-amplitude 
limit to the nonlinear center (0, u_). For v- in the range [0.55, 1) or (—1, —0.55], 
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there appears a single homoclinic orbit with vertex at v_ , hkewise enclosing a one- 
parameter family of periodic orbits about a nonlinear center {vc,U-). Together, 
these constitute all periodic solutions of (7.4) with s = for the double-well case 
(7.2). 

We may deduce by energy considerations that traveling periodic solutions of 
(7.1) must in fact be stationary, so that the Hamiltonian solutions just described 
constitute all periodic solutions of the van der Waals system (7.1). For, as noted 
in [Se.1-2], the flow of (7.1), for periodic solutions of period T, serves to decrease 
the mechanical energy 

(7.7) E{v):^E{v)+ I \u'^ dx = E{v) + ! \{u + eiv^fdx, 

Jo ^ Jo 

where 

(7.8) E{v):^J (^^\v'f + P{v)^dx 

denotes the Cahn Hilliard/van der Waals energy for the associated equilibrium 
problem. More precisely, 

(7.9) dE/dt = -{ei+£2) \ux\'^dx, 

Jo 

as may easily be derived from (7.3). Thus, we may deduce that u = constant for 
any periodic traveling-wave solution of (7.1), from which we immediately obtain 
= Vt = —sVx- It follows that either v is identically constant, in which case u is 
as well, or else s = 0; in either case, we find as asserted that the wave must be 
stationary. 

Moreover, one readily finds for = that any solution of (7.3) must satisfy the 
Euler-Lagrange equations 

(7.10) eie2Va:x = p{v-) - p{v) 
for the equilibrium variational problem 

(7.11) inmE{v), 

V 

yielding back the stationary-wave ordinary diff'erential equation (7.4) with s — 0. 

Stability. Using the result of Corollary 6.5, wc readily obtain: 

Theorem 7.1. Periodic orbits of (7.1) are unstable whenever dT /da > 0, in par- 
ticular, in the large- amplitude limit as they approach either a bounding homo- 
clinic or two-cycle. Orbits for which dT/da < are unstable if the spectrum of 
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(dum/du-) ^df{u-) is nonreal. In the small- amplitude limit as orbits approach a 
nonlinear center, they are unstable regardless of the value of sgn dT/da. 

Proof. Clearly, (H0)-(H2), hold, whence the results of Corollary 6.5 apply 
whenever (H3') holds, or dT/da ^ 0. By the final assertion of Corollary 6.5, we 
may without loss of generality take det{dum/ du^) ^ 0. The second assertion then 
follows by either of (6.28), (6.29), the former because in this case the determinant 
of the real matrix (dum/du-)~^df{u-) must be positive, as the product of com- 
plex conjugates. To establish the first assertion, we make use of the structure of 
equations (7.1). 

By direct computation, we obtain 

(7.12) =(^,(»_j -1). 

Moreover, from translational invariance of (7.1) with respect to w, in combination 
with the reflective symmetry u —u, x — > —x, we find without any computation 
that dum/du- must be of diagonal form 

(T.13) (« ;). 

with a single undetermined quantity a 7^ 0. Thus, assuming that AeiipUm/ du-) ^ 
0, or a 7^ 0, we find that the matrix 

(7.14) {dum/du.)-'df{v.,u.)=(^^,^^_^ -^q/") 

is trace-free, whence we may conclude that its spectrum is real if and only if its 
determinant is negative. But, this implies that (6.28) and (6.29) are mutually 
exclusive when dT/da > 0, yielding spectral instability. 

Finally, the third assertion follows simply by continuity of the Evans function, 
without invoking (H3') or Corollary 6.5, once we observe that the limiting, constant- 
coefficient equations at a nonlinear center are always unstable, since the trace-free 
real matrix df{vc,Uc) must necessarily have pure imaginary eigenvalues (or else 
(vc,ttc) would instead be a saddle). □ 

Remark 7.2. In the special case of the double-well potential (7.2), monotonicity 
dT/da > holds for all periodic van der Waals orbits. This follows as a straight- 
forward application of a monotonicity theorem of Schaaf (see [Sc, p. 102], or [LP, 
Theorem 7.1, p. 331]), which asserts that monotonicity holds for the nonlinear 
oscillator 



(7.15) 



v" + iJi{v) = 0, 
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provided: (i) jj,' > implies 5(/i")^ — Sn'/j,"' > 0; and, (ii) n' — implies n/j," < 
0. Noting that the van der Waals traveling- wave ordinary differential equation 
is of form (7.15), with fi{v) = p{v) — p{v-), we readily verify (i) and (ii) in the 
double-well case p{v) = 2v{l — f ^), by direct computation. In particular, note that 
(i) follows simply from the fact that p"'{v) — —12v < 0. This finding is quite 
significant, since it shows that there is at least one model exhibiting oscillatory 
pattern formation for which no stable periodic waves exist, thus indicating the 
presence of some alternative mechanism for the formation of oscillatory patterns. 

Remark 7.3. By Remark 5.10, in conjunction with Lemma 6.4, we find that 
small-amplitude periodic waves are stable with respect to periodic perturbations (of 
the same period T) if and only if dT/da > 0; in particular, they are stable for the 
double- well potential, as discussed just above. Likewise, numerical evaluation of 
dum/du- for large-amplitude waves in the two-cycle configuration yields F > 0, 
suggestive of stability under T-periodic perturbation, a conjecture that is further 
supported by numerical approximation of the {k = 0)-spcctrum of L; see Appendix 
A. On the other hand, we have already seen analytically that large-amplitude waves 
in the homoclinic case are unstable with respect to T-periodic perturbations. 

These results are quite interesting from the point of view of the associated vari- 
ational problem (7.11). For, they imply that periodic solutions may or may not be 
stable critical points of the energy E, depending on the details of their structure. 
This is in sharp contrast to the results of [CGS.1-2,GM] for the same problem 
with Neumann in place of periodic boundary conditions, in which critical points 
are seen to be stable if and only if they are monotone. In both cases, stability is 
interpreted with respect to perturbations conserving mass. Note, in the conserva- 
tion law setting, that stability is at best orbital, in the sense of convergence to the 
three-parameter family of nearby periodic solutions of the same period T, with the 
two conserved quantities given by the mass of the perturbed solution determining 
the limiting solution up to phase shift (translation). 

To put things slightly differently, stability in the periodic setting is not deter- 
minable as in the Neumann setting by Sturm-Liouville-type considerations alone, 
but must indeed be estimated by some such computation as we have carried out. 
Similarly, stability for the variational problem on the whole line is more subtle than 
that for the (Neumann) problem on the bounded interval, requiring further analysis 
as in [Z.l]. 

Quadratic-flux models. Planar systems (1.1) with quadratic flux / have been 
studied as qualitative models for multiphase flow in porous media near an "ellip- 
tic boundary" where characteristics of the flux Jacobian df{u) coalesce; see, e.g., 
[SSh,MP]. Here, following [GZ], we consider the restricted class of 2 x 2 systems 
with / quadratic and B constant for which, additionally, the traveling-wave ordi- 
nary differential equation 

u'^Biu)-\fiu)-fiu-)-s{u-U-)] 
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is Hamiltonian, i.e., 

tr{B-\f'{u)-sI)) = 

for some choice of s. This is a codimension-2 subclass of the quadratic-flux models, 
as parametrized by the coefficients of B and /. 

Generically, such systems can be reduced by an affine change of variables to the 
canonical form 

(7.16) + ^sh' 

where e — ±1. For e = 1 and s = 0, the traveling- wave ordinary differential 
equation then reduces to the Hamiltonian system 



f 

(7-17) , , , 

\V j \ UV I \ U-V 

with Hamiltonian 

(7.18) H{u,v) = -^ii^v^ — u^v + v'^) — {-v'^v — -u^v + V-V — U-V-u) = H{uq,vq) 

preserved along orbits [GZ]. 

Explicit computation of the phase portrait of (7.17) using (7.18) reveals the 
existence of a unique three-cycle configuration for the special parameter values 

(7.19) («_, ^_) = (±V3/2, 0), (0, -1.5), (0, -.5), 

corresponding, respectively to the three saddle equilibria at the vertices of the three- 
cycle, and the single nonlinear center contained within. Enclosed within the three- 
cycle is a nested one-parameter family of periodic orbits. The heteroclinic orbits 
making up the three-cycle lie on straight lines, forming an equilateral triangle; in 
fact the entire phase portrait has the same triple symmetry about the center as 
does the bounding triangle. The three states (7.19) at the vertices of the triangle 
may be regarded as analogous to 'pure phases in physical models for three-phase 
flow, and the interior of the triangle to the physical state space of volume fractions 
thereof [AMPZ.3]. 

This configuration bifurcates as (u_,w_) is varied into various one- and two- 
cycle configurations, each cycle likewise enclosing a nested one-parameter family of 
periodic orbits. Where dT/da ^ 0, these are the only nearby periodic orbits to be 
found in parameter space, where s is now allowed to vary; however, we have not so 
far been able to rule out the possibility of further periodic orbits with large speed 
s. 

In numerical experiments carried out in [AMPZ.5], systems (7.16) were seen to 
exhibit the same sort of oscillatory pattern formation discussed in [FL.1-2,CP], with 
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patterns composed of motifs with (w, v) profiles lying near either the three-cycle, 
or the various two-, and one-cycle configurations that bifurcate from it as {u-,V-) 
is varied from (7.19). Thus, the relevant issue from the point of view of pattern 
formation seems to be the stability or instability of large-amplitude periodic waves 
for parameters near the three-cycle values (7.19). Regarding this problem, we have 
a complete solution; indeed, the additional symmetries of the three-cycle case allow 
us to obtain results for general amplitude waves: 

Theorem 7.4. For B constant, symmetric, and positive-definite, periodic orbits 
of (7.16) are unstable for {u-,V-) lying near the three-cycle values (7.19), inde- 
pendent of amplitude, provided that dT/da ^ 0. Likewise, so long as dT /da ^ 0, 
large- amplitude waves are unstable for {u-,V-) corresponding to a homoclinic con- 
figuration in the phase portrait of (7.17); moreover, small- amplitude waves are un- 
stable for any phase configuration. 

Proof. Clearly, (H0)-(H2), hold, whence, again, the results of Corollary 6.5 
apply whenever (H3') holds, or dT/da ^ 0. Taking (tt_,z;_) to be the nonlinear 
center (fXc Vc) = (0, —0.5), we find using the triple symmetry of the phase portrait 
that d{um,Vm)/d{u-,V-) must be a multiple of the identity, for any periodic orbit 
whatsoever. More precisely, we note that, for {u-,V-) = {uc,Vc) = (0,-0.5), 
traveling-wave ordinary differential equation (7.17) has reflective symmetry about 
each of the lines from {uc, f c) through the vertices of the surrounding three-cycle; 
moreover, each individual symmetry is preserved as {u-,V-) is varied along the 
respective line of symmetry. It follows that, under perturbations of {u-,V-) in 
direction u, where u is the direction vector of a line of symmetry, must be confined 
to the same line of symmetry, i.e., the resulting perturbation of the mean {um, Vm) 
must lie in the same direction for any choice of initial periodic orbit (w, v). In 
other words, each of the direction vectors of the three lines of symmetry is an 
eigenvector of d{um,Vm)/d{u-,V-), whence d{um,Vm)/d{u-,V-) must be a (real) 
multiple of the identity. 



Thus, 



df{u-,v-)d{um, Vm)/d{U-,V-) 



-1 



is simply a real multiple of 



(7.20) 




But, this matrix is similar to the anti-symmetric matrix 




hence must have pure imaginary spectrum. (Note: in particular, it is again trace- 
free, as in the proof of Theorem 7.1.) We thus obtain instability by (6.29), for 
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waves of arbitrary amplitude. The result for V-) near (0, —0.5) follows by the 
same argument, together with continuity of d{um,Vm)/d{u-,V-). 

The large-amplitude homoclinic result has been established already by several 
different techniques. Likewise, the small-amplitude result follows from the instabil- 
ity of the constant coefficient limit, as discussed previously in Remark 4.5. □ 

Remark 7.5. The proof of Theorem 7.4 shows also that, for (u-,V-) near the 
3-cycle values (7.19) and dT/da > 0: - in particular, for large-amplitude periodic 
orbits approaching the 3-cycle, the stability index F is positive, suggestive of sta- 
bility with respect to T-periodic perturbations. This conjecture is supported by 
numerical approximation of the (k = 0)-spectrum, which also indicates stability 
with respect to T-periodic perturbations of large-amplitude periodic solutions in 
the 2-cycle limit; see Appendix A. This is in contrast to the T-periodic instabil- 
ity observed in the large- amplitude homoclinic limit for both van der Waals and 
quadratic flux models. 

Remark 7.6. The general features of the global phase portrait structure de- 
scribed above for the Hamiltonian models considered here have been shown in 
[AMPZ.4] to hold, more generally, for all 2 x 2 quadratic-flux models of type I un- 
der the viscous classification scheme of Hurley [HP] . In particular, there exist unique 
{u-,V-, s) such that the phase portrait of the traveling- wave ordinary differential 
equation contains a three-cycle, of which the connecting heteroclinic orbits must lie 
along straight lines; moreover, it can be shown by affine reduction to canonical form 
that the corresponding ordinary differential equation is quasi-Hamiltonian [CL]. It 
follows that the three-cycle encloses a nested one-parameter family of periodic or- 
bits, as in the "pure" Hamiltonian case considered above. Likewise, formation of 
oscillatory patterns has been observed numerically for such systems [AMPZ.5], with 
patterns again lying near the special three-cycle configuration. 



Appendix A. Numerical experiments 

In this appendix, we describe numerical methods for the location of the zero 
set of the periodic Evans function D{k, A), or equivalently the spectrum of a linear 
operator L with periodic coefficients, and, in particular, the well-conditioned nu- 
merical verification of instability. We demonstrate these methods for the example 
systems considered in Section 7, amplifying and (in the case of the quadratic flux 
models) extending the analytical results there obtained. The real advantage of the 
numerical approach, of course, is that it is applicable to general models. 

The algorithm. We use the following basic strategy for numerical verification 
of instability: 

1. Determine the range of existence of periodic orbits, and (numerically) solve 
for a representative sampling of profiles. 
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For each profile: 

2. Find spectral curves A(/c), i.e., zero-level sets of -D(A, /c); 

3. Choose k = k* giving a maximum (positive) value of Re {X{k)) and a numer- 
ically advantageous closed contour F contained in the strictly positive half-plane 
Re A > and containing A(/c*). 

4. Perform a winding number calculation in A, with k — k* held fixed, to 
establish that there indeed lies a root of D{k*,-) in the vicinity of the original 
numerical approximation of \{k*) in steps 2 and 3. 

The point of steps 3 and 4 is that the curves found in step 2 are "graphical" 
approximations of the zero-level sets, which, though extremely illuminating, do not 
give conclusive information about existence of zeroes. Step 4 by contrast is both 
accurate and well-conditioned, and could be used as the basis for numerical proof; 
however, it requires an a prior guess for A;*, X{k*). 

Choosing k* and F. By a standard Carding type energy estimate on (2.6), the 
zeros of the Evans function must lie in the truncated wedge 

y - {A e C : Im (A) + Re (A) < r, Im (A) - Re (A) > -r. Re (A) < r/4} 

where r = ||/(0lloo/ll-^ll2) /(C) = ll^lOlb- Obviously, this is not a bounded region 
but it does give us a bound on the unstable zeros of D. We find spectral curves/zero 
sets of D numerically and observe that those are contained in V. We can see also 
agreement between our numerics and the analysis in Sections 1-7 (see discussion 
under Applications, below). 

It is more tricky to find the zero sets of the two-parameter periodic-type Evans 
function than of one-parameter front-type Evans function. Adding to the difficulty 
is the fact that -D(A, k) vanishes to {n + l)st order in A at (0,0), where n is the 
dimension of variable u, so that there are n + 1 spectral curves bifurcating through 
the origin. The apparently straightforward step 2 in the algorithm is thus in fact a 
tricky problem in graphical display. To find spectral curves/zero sets of D, we use 
Miiller's method which can approximate complex roots and converges to the root 
for any initial approximation choice. We also use defiation to achieve a reasonable 
result near the bifurcation point (/c. A) = (0,0). 

Since the argument principle will be used to count the zeros inside some contour 
F, care must be taken that F itself does not contain any zeros. We avoid the origin 
and high multiplicities by looking at A; 7^ 0. A natural choice is to choose the 
value of k which gives the maximum value of the positive Re (A) from the spectral 
curves, and for simplicity choose F as the rectangle whose center is that maximum 
positive value of Re (A) in the right half plane. In this process, we find that a 
practical upper bound on the size of the contour is enforced by the need to obtain 
good error estimates, which degrade rapidly for high frequencies. This particular k 
allows good numerics. When the computed winding number exceeds zero, the wave 
being tested is unstable; otherwise the wave is stable. 
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Applying the Argument Principle. Here we follow the approach of [Br. 1-2]. The 
first problem in calculating the winding number of D[T] lies in evaluating D{X, k) for 
a fixed value of k. We use the maximal value k = k* described above. The second 
problem lies in partitioning F into subcontours such that -D[Fi] lies in a slit 
plane; a slit plane being a subset of C of the form C\{ra : a E C, < r E M.} . 
This will ensure that the calculated winding number is correct. We partition F so 
that -D[Fj] is contained in a half plane whose boundary passes through the origin. 

For w,z G C, define {w, z) to be the dot product between the real 2-vectors 
(Re (t(;),Im (w)) and (Re (^),Im (z)) : 

(8.1) {w, z) = Re (M;)Re {z) + Im (M;)Im {z). 

If T{si) is the starting point of Fj, it suffices to find s^+i and (3 e [0, 1) such that 

(8.2) {D{s, + ^s),D{si))> P\D{si)\ 

for < As < Si_|_i — Si- In other words, the component of D{si + As) in the 
direction of D{si) should be a positive multiple of D{si). Strictly speaking, the 
right side of inequality (8.2) could simply be in practice. However, to allow for 
numerical error, a positive value of /? should be chosen. We choose (3 = 1/2 over 
all. This condition allows checking whether a given partition is appropriate. 

Error Estimates. In contrast to the front case [Br. 1-2], computation of D for 
moderate values of |A|, given exact values for coefiicients of L, is straightforward and 
numerically well-conditioned, and for practical purposes may be regarded as exact. 
The main approximation error comes rather from the computation of coefficients A 
and B in (1.4) through solution of the traveling- wave ordinary differential equation 
(1.8), which is likewise a standard and relatively well-conditioned problem; here we 
use a fourth order Runge-Kutta method, for which analytical error estimates may 
be foimd, for example, in [SB]. For large |A|, it is rather the 0(|A|) exponential 
growth rate in the eigenvalue ordinary differential equation that dominates the 
growth of errors; indeed, it is this consideration that effectively limits the size of 
our contours in the winding number calculations described above. However, again, 
the relevant analytic error bounds are standard and straightforward. 

We therefore omit the discussion of analytical error estimates, instead carrying 
out a numerical convergence study. Specifically, we use Richardson extrapolation 
to compute "exact values" to which we compare when we calculate relative errors 
and to check the convergence of D as the step size h becomes small. The resulting 
estimated errors are seen to be quite small, consistent with the above discussion. 

Numerical experiments. We now describe the results of our numerical exper- 
iments for the two classes of systems discussed in Section 7. 

Van der Waals gas or solid mechanics. Our first set of experiments is for the 
viscous-capillary p-system (7.1) with standard double-well potential (7.2). As dis- 
cussed in Section 7, periodic solutions of this model are of necessity stationary. 
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satisfying the Hamiltonian system (7.4) with s = 0, and this system features for 
v_ = ±1 a one parameter family of heterochnic two-cycles, indexed by U- due to 
Galilean invariance in the speed u, in which there lie nested families of periodic 
orbits. Likewise, there exists a two parameter family of homoclinic cycles for v- 
in range [0.55, 1), and U- again arbitrary by Galilean invariance, in which there 
lie nested families of periodic orbits. Without loss of generality (using transla- 
tion invariance), we may fix tt_ =0, reducing our study to a single 2-cycle and a 
one-parameter family of homoclinic orbits, and their enclosed periodic orbits, see 
Figures 6A-B below. 

First, we use a fourth order Runge-Kutta method to get a profile {v, u) by solv- 
ing the Hamiltonian system. By plotting v versus tt, we can see clearly there exist 
periodic orbits (Fig. 2A) inside the two cycle. From these orbits, we observe that a 
monotone relationship between the periods and the amplitudes (Fig. 2B). Second, 
we use the method again to solve (2.7). Finally, we use the LU decomposition 
method to get a value of D. We also find periodic orbits and the monotone rela- 
tionship inside the homoclinic orbit numerically. We study several periodic orbits 
in each cycle. 



The orbits of Van der Waals eq. 



TTie period as a function of amplitude for Van der Waals eq. 





Figure 2a. Existence of periodic 
orbits. 



Figure 2b. Period as a function 
of amplitude. 



Figure 3 shows spectral curves of L, or zero sets of -D, for various periodic orbits 
within the 2-cycle depicted in figure 2A, starting with the (zero amplitude) constant 
coefficient case and shows how the spectral curves are changed as the amplitude 
becomes large. We see what differences are made by change of amplitudes. These 
figures are extremely illuminating, showing in a way that calculations can not the 
detailed evolution of the spectrum as the amplitude of the underlying periodic orbit 
is varied from zero (constant coefficient case) to the limiting amplitude as orbits 
approach their bounding polycycle. They both verify and amplify the observations 
obtained earlier by Evans function calculations. In particular, we immediately see 
that instability is rather dramatic, whereas the analytical results via Taylor expan- 
sion about the origin imply only that there exist infinitesimally unstable spectra; 
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The e-values of Van der Waals eq. -const, case 



e-values of Van der Waals eq. - non-consL case(amp=.066) 



.'oe 
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Figure 3. Spectral curves of the Evans function D. 



moreover, this instabihty is not confined to smaU or large amplitude limits, but 
persists across all amplitudes. 

This suggests the strategy described above for numerical verification of insta- 
bility; namely, to fix k at the value k* corresponding to the maximally unstable 
spectrum A* = A(/c*) obtained by graphical examination, then perform a winding 
number calculation on a well-chosen contour about A*. 

As for stability results. Figures 4 and 5 show F and D[r] for periodic waves 
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inside a two cycle and a homoclinic orbit respectively. Here, D[r] is normalized by 
the mapping D i — )• 0.1 * {i + 1) * D/\D\ for z = 0, 1, ...,39 (Here, i is the index 
of A on the contour F. The winding numbers are greater than 1 both cases so the 
periodic waves are unstable. 
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Figure 4. The number of unstable eigenvalues in the two-cycle case. 
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Figure 5. The number of unstable eigenvalues in the homoclinic case. 
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We next perform the error analysis that was mentioned earher. There is no 
analytic solution of the eigenvalue ordinary differential equation (2.7) for the van 
der Waals equations. We have to rely entirely on numerical estimates to gauge the 
accuracy of the Evans function method. In the tables below, third order polynomial 
extrapolation was used to compute the values for r = where r is a local truncation 
error depending on the step size /i, i.e., r := C/i^ for some positive constant C (recall 
that the error to lowest order comes from solution of the profile equation by the 
fourth order Runge-Kutta scheme). 



T 




i:>(0.32) 


Extrapolated Relative error 


10" 


■5 


3.358392 


0.0121955 


10" 


■6 


3.358958 


0.0120290 


10" 


■7 


3.359097 


0.0119882 


10- 


•8 


3.399841 


0.0000041 







3.399855 







Table 1. D{{).Z2) for the 


van der Waals equations 


T 


i:>(0.52 - O.H) 


Extrapolated Relative error 


10" 


^ -4.636345 - 12.13182z 


0.0082597282 


10" 


^ -4.636959 - 12.13066i 


0.0082772553 


10" 


^ -4.637134 - 12.13046i 


0.0006327496 


10- 


^ -4.722007 - 12.19752^ 


0.0000185195 





-4.721931 - 12.19775i 





Table 2. D(0.52 — O.lz) for the van der Waals equations 

Starting with the extrapolation results, examine the above Tables which describe 
an illustrative calculation when the amplitude of periodic orbit is 0.32. The values 
chosen for A are the points on V closest to and furthest from the center of F; in 
practice, these are found to give the extrema of errors. The extrapolated values are 
used as exact values in computing extrapolated relative errors. Note that r = 10~^ 
is sufficient for calculating -D(A, k) with approximate estimated accuracy 0.05. 
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Similar calculations were carried out for all mesh points A on the contours for 
each periodic orbit under consideration. For example, in the two-cycle case, we 
tested all periodic orbits indicated in Figure 6A. The points in the Figure represent 
the amplitudes of each periodic orbit. Note: translation invariance greatly reduces 
the dimension of the problem. Likewise we consider homoclinic profiles with vertices 
V- e [0.55, 1) with vertex mesh 0.05 (see Fig. 6B). For each homoclinic case, we then 
computed periodics with amplitude mesh 0(0.01) where amplitude e [0,0(0.1)]. 
In each case, the winding number was found to be greater than 1. The bounds for 
the winding number calculations are found conservatively to be order 10~^, i.e., the 
worst case Rouche's error over all calculation was 



max- 



|AD| 



< O10"2 « 1/2, 



for some constant C (Table). It is comfortably below the bound 1 needed to apply 
Rouches Theorem. Thus, we have numerical instability across the whole parameter 
range of existence. 





Figure 6a. The plotted points rep- 
resent the amplitudes for which the 
experiments were carried out in the 
two-cycle case. 



Figure 6b. The plotted points rep- 
resent the amplitudes for which the 
experiments were carried out in the 
homoclinic case. 



As discussed in the main body of the paper, there are various interesting tran- 
sitions of the number of unstable 0-eigenvalues of the periodic waves with respect 
to different amplitude. Recall that the small amplitude limit is the constant co- 
efl&cient case. For fixed /c = 0, the Evans function has four roots at the origin if 
the amplitude is zero, and three roots at zero otherwise when n = 2; see Remark 
4.5. But, we can track a fourth small root of D as the amplitude changes. When 
the amplitude is small, the fourth root is typically in the left half plane of Re (A) 
(Remark 4.5). In the homoclinic the amplitude becomes large, the fourth 
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-0.02 -0.01 0.01 0.02 -0.1 -0.05 0.05 0.1 



lambda lambda 

Figure 7. Transition of the number of unstable 0-eigenvalues for the 
van der Waals equations, homochnic case. 



root moves from the left to the right half plane crossing the origin (Fig. 7); see 
Remark 5.11. For the two-cycle case, no such transition occurs. 



The quadratic-flux models We next consider the class of quadratic flux models 
(7.16) with identity viscosity matrix B = I, applying the same methods that we 
used for the van der Waals equations. First, we get a profile {u, v) by solving (7.17). 
By plotting u versus v, we can see clearly that there exist periodic orbits (Fig. 8A). 
From these orbits, we observe that there is a monotone relationship between the 
periods and the amplitudes (Fig. 8B). Finally, we use the LU decomposition method 
to get a value of D. We have spectral curves of A with respect to various amplitudes 
(Fig. 9). 

The quadratic-flux model admits, for s = 0, a single 3-cycle solution, a 1- 
parameter family of 2-cycle solutions, and a 2-parameter family of homoclinic so- 
lutions, each containing a 1-parameter family of nested periodic orbits. We restrict 
to the case s near zero, for which these are the only periodic solutions. The lo- 
cation of 3- and 2-cycles may be determined analytically; however, a new issue as 
compared to the van der Waals case is that the region of existence for homoclinics 
is no longer analytically obtainable. Note: The regions of existence for two-cycles 
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The orbits of the quadratic flux model Tlie period as a function of amplitude tor quadratic flux model 




Figure 8a. Existence of periodic Figure 8b. Period as a function 

orbits. of amplitude. 




Figure 9. The spectral curves of Evans function D. 



and three-cycles are analytically obtainable, as discussed in [AMPZ.3] and [GZ]. 
We determine the region of existence of homo clinic orbits numerically by varying 
{u-,V-) where U- = (—0.86, 0.86), f_ G [—0.5,0.5). Figure 10 shows the break-up 
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Figure 11a. The plotted points 
represent the vertices of one- and 
two-cycle for which the experiments 
were carried out. 



Figure 11b. The region where ver- 
tices of homoclinics lie. 



of the homoclinic orbit of this model. Figure 11 A shows the tested points 
which are vertices of homoclinics (•) and two-cycles (*), respectively. The vertices 
of the (unique) three-cycle are the three vertices of the surrounding triangle, the 
vertices of two-cycles lie on the edges, and the vertices of homoclinics lie within 
shaded region (Fig. IIB). For homoclinics, we only need to treat the region of 
{u-,V-) inside the triangle since every outside point has a corresponding inside 
point yielding the same dynamical system (i.e., stationary points occur in pairs, 
one inside, one outside). 

Inside each polycycle, we find several periodic orbits. To these, we apply the 
error analysis and the winding number calculations as for the previous model. The 
same error analysis applies here as in the previous section so we shall not discuss this 
aspect. The maximum Rouche's error over the entire region of testing was C10~^ 
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for some positive constant C, where in each case the winding number is greater 
than 1. Figures 12, 13, and 14 show the instabihty of each representative periodic 
orbit within thrcc-cyclc, two-cycle and homochnic. We see that small amplitude 
orbits and also large amplitude orbits lying near the three-cycle are unstable, as 
shown in Theorem 7.4, and also intermediate amplitude orbits are unstable. Thus, 
for this model system, too, we have instability across the whole region of existence. 

Finally, similarly as in the case of the van der Waals equations, we may track the 
number of unstable 0-eigenvalues of the periodic waves with respect to amplitude 
in the homoclinic, two-, and three-cycle cases. The result in the homoclinic case 
is that there is a transition from zero to one unstable root as amplitude increases. 
In the two- and three-cycle cases, there is no such transition, indicating that waves 
of all amplitudes are stable under T-periodic perturbations, even though they are 
unstable under perturbations of general period. (Recall, T-periodic stability can 
be shown analytically for small- amplitude waves by consideration of the constant- 
coefficient case; see Remark 4.5.) 
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Figure 12. The number of unstable eigenvalues, three-cycle case. 




Discussion 

The instability of stationary periodic solutions of the example systems of Section 
7 has been shown numerically. The code is general and easy to manipulate. It 
would be quite interesting for more general systems to either establish analytically 
a corresponding result of uniform instability or find a counter-example of a stable 
waves. Note that the code can be used also to investigate stability (for all k, not 
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Gamma Contour The number of e-values (3), amp=.2 
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Figure 13. The number of unstable eigenvalues, two-cycle case. 
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Figure 14. The number of unstable eigenvalues, homoclinic case. 

one k), though it has some difficulties near k — due to numerical sensitivity 
associated with the multiple root in A occurring there. This allows the systematic 
exploration of periodic wave stability for families of conservation laws, and also for 
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Cahn Hilliard and reaction diffusion equations. 



Appendix B. Stationary phase approximation 

In this appendix, we reproduce some unpubhshed computations of [AMPZ.5] 
referred to in the introduction, concerning the linear response under perturbation 
of a linearly unstable constant solution tt = tto of a conservation law with viscosity 
(1.1). Consider the linearized equation (1.4), where A and B are now constant. 
Taking the Fourier transform, we obtain 

(10.1) vt = P{ik)v := {-ikA - k^B)v. 

The fundamental solution of (10.1) is clearly ^'(/c,t) = QP{ik)t^ whence the funda- 
mental solution of (1.4), i.e., the solution with initial data given by a Dirac delta 
function centered at the origin, is given by the inverse Fourier transform 

(10.2) ^{x,t) = — e'^'^e^^'^'^Uk, 

and the Green function for (1.4) by G(x, t\y) — "^{x — y, t). 

The n eigenvalues of the symbol P{ik) determine n dispersion relations 

(10.3) A = A,(A;), j = l,...,n, 

with linearized stability corresponding to the condition Re < for all /c G M. 
Denote by k* the critical frequency at which Re A-, takes on its maximum value. 
For simplicity, we restrict to the generic situation that Xj is isolated from the rest of 
the spectrum of P{ik) at k = kj, and Xj{k*) is a nondegenerate maximum. Thus, 
Xj admits a second-order Taylor expansion 

(10.4) Xj{k) = Xj{k*) - iaj{k - k*) - (3j{k - k*f + . . . 

about /c*, where aj is real, and Re 13 j > 0. Let Rj{k) = Rj{k^) + . . . and Lj{k) = 
Lj {k^) + . . . denote right and left eigenvectors and Ilj := RjL* / {Rj, Lj) the spectral 
projection of P associated with Xj. 

Decomposing (10.2) by spectral resolution as '^{x,t) = J2j where 

1 

(10.5) '^j{x,t) ■=Y1^J e''"=e^^^''^'Ujik) dk, 

we may estimate each of the by stationary phase approximation about kj as 
(10.6) 

f+OO 



-1 I— too 

*j(a;, t) ~ := — / e*'=^e(^^-('^P-^"^-('=-'=P-^^('=-'=^*)')%(A;*) dk 

27'' J — oo 
1 r+oo 

27r J— oo 
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Let US now restrict to the simplest case B = I. In this situation, the dispersion 
relations are exactly Xj{k) = —ikaj — k^, where aj are the eigenvalues of A, and 
the associated eigenprojections IT-, are just the eigenprojections rjlj/{rj,lj) of A 
associated with aj, where rj and Ij denote associated right and left eigenvectors of 
A. Thus, linear stability is equivalent to reality of aj, or hyperbolicity of the first- 
order coefficient A. By direct calculation, we find that kj = Im aj/2, aj — Re aj, 
X{kj) = (Im — zRe a^Im aj/2, and so 

(10.7) aj = -Im \{k*)/k*. 

Likewise, noting that I3j is real, we may explicitly carry out the Fourier inversion 
in the final factor of (10.6) to obtain 



X 



(10.8) ^gReA(fc;)f 'J'j 



rA* Jk*{x-ajt) -{x-ajtf/4t 



X 



{rj,lj) ViTrt 



= g(Im aj2yt '3" 3 ^ 



2 r l* g«(Im aj/2)(a;— Re ajt) ^—(x—Ke ajt) /4t 



where we have used (10.7) in the second equality. That is, the behavior of $ in 
the jth mode is, to lowest order, that of a time-exponentially growing modulated 
Gaussian wave packet traveling with speed aj, with exponential growth rate equal 
to the maximum real part Re Xj{kj) of Xj, and modulating spatial frequency equal 
to the critical frequency kj. 

In the case of a linearly stable mode Im aj = 0, (10.8) reduces to the usual 

Gaussian "linear diffusion wave" -^^r-fry x ^/j^^ of [LZeJ. By analogy, we 

might call (10.8) an "unstable linear diffusion wave." Note that unstable modes 
Im aj 7^ group in conjugate pairs, with equal propagation speeds aj = Re aj, and 
thus their diffusion waves combine to give a single real, sinusoidal diffusion wave 
of dimension two. The proper analogy is, therefore, to the "generalized diffusion 
wave" defined in [LZe] in the case of repeated real eigenvalues of A (the boundary 
case between strict hyperbolicity and ellipticity of A). 

Example 10.1. For the equations of van der Waal gas dynamics or elasticity, we 

have B = I and A = df{uo, ~ (^p'^y ^ 0^ ) ' eigenvalues aj = ±y/—p'{vo). 

Thus, linearized stability of the constant solution (u,v) = {uo,vo) is equivalent to 
p'{vo) < 0, and failure of stability entails a pair of pure imaginary eigenvalues 
Oj = ±ir. In this case, therefore, there is a single conjugate pair of standing 
unstable linear diffusion waves, aj = 0, with time-exponential growth rate e'^ 
and modulating spatial frequency e*"^^/^. 
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Application to pattern formation. Now, let us consider the implications of 
the linear estimate (10.8) on the behavior under perturbation of solutions of the 
nonlinear equations, for initial data approaching a linearly unstable constant state 
Uoo as X — >^ +00. In view of the global nonlinear structure of our models, let us 
make the crude approximation that nonlinear effects take over at some prescribed 
amplitude A, turning off the time-exponential growth and converting the sinusoidal 
oscillations of the linearized approximation (10.8) to a series of traveling fronts mov- 
ing with approximately zero speed (the last assumption implies some normalization 
to achieve a rest frame). Looking at a single unstable mode j, the picture that 
emerges from this crude model on the half line [0, +00) is an oscillatory pattern 
of phase transitions separated from the constant state by a modulated front whose 
location x{t) is determined by the property that 

(10.9) \^{xj{t),t)\^A, 

i.e., ignoring the algebraic factor t~^/'^ in (10.8), travels with approximate speed 
aj determined by (Im aj/2)^t — {ajt — ajt)'^/4t = 0. Solving, and recalling that 
aj = Re aj, Re X{kj) = (Im aj/2)^, we obtain a value of 

(10.10) CTj = Re ttj + |Im aj\ 

for the approximate front speed. 

Since waves under this model are "born" in the linearized regime, we can estimate 
the approximate wave-length inherited in the nonlinear regime by dividing the 
distance d = ajt swept out by the front x{t) from time to time t, divided by 
the number of peaks of the linearized unstable diffusion wave (10.8) that it has 
crossed during that time. Working in the rest frame x = ajt of the diffusion 
wave, we find that the number of peaks is {aj — aj)t divided by the wavelength 
2n/k* — 27r/Im aj/2 of its spatial modulation, yielding an estimated wavelength of 



(27r/fe*) 



This predicts an average (Doppler-shifted) frequency for the nonlinear pattern of 



aj - aj - 


) 


aj , 




Im 




Re aj -\- 


Im a J 



(10.11) 



In the (standing) case a-, = Re aj = of Example 10.1, this is just the critical 
frequency k* = Im aj/2. Note: here and above we are implicitly assuming that the 
front speed aj is positive, i.e., the pattern is expansive. This is in fact a necessary 
condition for pattern formation to occur. 
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Numerical experiments in [AMPZ.4] showed close agreement between the pre- 
dictions of (10.10) and (10.11) and behavior of actual patterns in solutions of (1.1). 
Likewise, the condition aj > was seen to well predict nonlinear instability, and the 
appearance of pattern formation, as model parameters (e.g., Riemann endstates) 
were varied. 

Remark 10.2. In the linearly stable case, Liu and Zeng [LZe] derive also a self- 
similar nonlinear diffusion wave refining the linear estimate (10.8), by judiciously 
appending to an approximately diagonalized version of the linearized equations 
(1.4) the "diagonal" quadratic-order terms in the Taylor expansion of the nonlinear 
flux /. This approximation makes essential use of the small-amplitude nature of 
solutions. In the large- amplitude linearly unstable case, we have no such simple 
description of a nonlinear diffusion wave; indeed, the relevant entity seems to be the 
entire metastable pattern described above. We might neatly summarize the conclu- 
sions of this paper (together with those of [AMPZ.4-5]) as suggesting that patterns 
observed in dynamical phase-transitional models represent metastable nonlinear 
diffusion waves rather than stable periodic states. 
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